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1.  Introduction 

During  the  period  of  April  1976  through  April  1978  DPSI 
was  under  contract  (F19628-76-C-0241)  to  the  Convective  Cloud 
Physics  Branch  (LYC) , Meteorology  Division  of  the  Air  Force 
Geophysics  Laboratory  (AFGL) . The  purpose  of  this  contract 
was  to  develop  and  apply  mathematical  procedures  to  a variety 
of  standard  and  non-standard  cloud  physics  research  data. 

The  work  performed  under  this  contract  has  been  submitted 
in  detail  in  a set  of  23  monthly  reports.  Rather  than  summari- 
zing these  here,  this  report  includes  a full  set  of  updated 
documentation  and  operating  instructions  for  the  programs  de- 
veloped under  this  contract. 
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2.1  KN1UTIL 

2.1.1  Program  1 :scription 

Program  KN1UTIL  was  written  to  examine  or  copy  a PMS-lD 
data  tape  to  verify  the  correct  operation  of  the  associated 
hardware.  It  also  provides  an  input  summary  that  is  invalu- 
able to  the  checkout  of  the  post  processing  programs  that 
utilize  this  data. 

KN1UTIL  is  usually  executed  with  cards  through  the 
CDC  6600  batch  processor  and  will  produce  either  an  output 
listing  or  a magnetic  tape.  It  can  list  the  input  data  in 
various  ways: 

1)  an  octal  dump  of  each  record 

2)  a decimal  dump  of  each  record 

3)  a decimal  dump  of  any  selected 
probe  or  VCO  string 

4)  a decimal  listing  of  the  PMS  status 
words 

Each  method  gives  a multi  file  reading  capability  and  also 
a means  to  skip  records  when  listing. 

An  additional  facility  of  KN1UTIL  is  a copy  option  that 
allows  various  records  from  different  files  on  the  PMS-lD 
tape  to  be  copied  to  a new  tape  in  the  same  format. 

The  input  data  tape  can  come  from  three  sources:  the 
Kennedy  recorders  on  the  C130-A  or  C130-E  aircraft  or  from 
the  Learjet  PMS-lD  device  after  being  preprocessed  by 
LEARPMS. 
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2.1.2  CONTROL  CARDS 

JOBID, CM60000 , T100 , TP2 . * Prob.  No.  Name 

VSN , TAPEl=TAPENO , TAPE2=TAPENO. 

ATTACH  LGO , KN1UTILBIN , ID=GLASS , MR= 1 . 

REQUEST, TAPE1 , S , HI . 

REQUEST , TAPE 2 , S , HI . 

FILE (TAPEl , RT=U , BT=K , MRL=102  4 , MBL=1 02  4 , RB=1 , BFS=105 ) 
FILE (TAPE2 , RT=U/BT=K,MRL=1024 ,MBL=1024 , RB=1 , BFS=105 ) 
LDSET ,FILES=TAPE1/TAPE2 , PRESET=ZERO . 

MAP , OFF . 

LGO. 

REWIND , BPARAM , CP ARAM , DP ARAM .* * 

COPY, BP ARAM. ** 

COPY , CPARAM . * * 

COPY , DPARAM . ** 

7/8/9 

DATA  DECK 

7/8/9 
6/7/ 8/9 


If  not  duplicating  a tape 

1.  change  TP2  to  TP1 

2.  remove  TAPE2=TAPENO  from  VSN  card 

3.  remove  REQUEST  TAPE 2 card 


**  Remove  the  cards  if  status  word  dump  is  not  desired 
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2.1.3  DATA  CAR  '3 

1.  Decimal  I imp  - one  data  card 

cc  i 1 DEC 

cc  5-7  number  of  files  dumped 

cc  8-10  record  indicator 

i.e.  1 = every  record 

2 = every  other  record 
6 = every  sixth  record 

ex.  to  dump  decimally,  two  files  every  10th  record 

-DEC — 2-10  - = blank 

2.  Octal  Dump  - one  data  card 

cc  2-4  OCT 

cc  5-7  same  as  decimal  dump 

cc  8-10 

ex.  to  dump  octally,  nine  files  every  100th  record 
-OCT — 9100  - = blank 

3.  Selective  Record  Dump  - n data  cards  (decimal  only) 
card  1 

cc  2-4  REC 

cc  5-7  (rj)  number  of  files  dumped 

cc  8-10  (rj)  record  indicator 

cards  2...n  (one  card  for  each  set  of  records) 
cc  1-6 (rj)  starting  record  to  dump 

cc  7-12 (rj)  ending  record  to  dump 

ex.  to  dump  every  10th  record 
from  125  to  350  and 
from  1000  to  1500 


-REC — 1-10 
125 350 

--1000  — 1500  - = blank 


4.  Tape  Copying  - n data  cards 
card  1 

cc  2-4  DUP 

cc  5-7  number  of  files  to  copy 

cards  2...n  (one  card  for  each  file  to  be  copied) 

cc  1-6 (rj)  number  of  records  to  skip  before  copying 

cc  7-12 (rj)  number  of  records  to  copy 

ex.  to  create  a tape  of: 

records  701  to  750  from  file  1, 
and  records  25  to  50  from  file  3, 
the  data  deck  is  ... 

-DUP — 3 

700 50 

0 0 

24 26  - = blank 

5.  Selective  Probe  Dump  (1  data  card) 


cc 

2-4 

DEC 

cc 

5-7 

number  of  files  dumped 

cc 

8-10 

record  indicator 

i.e.  1 = every  record 

2 = every  other  record 

6 = every  sixth  record 

cc  15  selected  probe 

BLANK/ 0 for  regular  decimal  dump  (OPTION  1) 

1 for  scatter  probe  only 

2 for  cloud  probe  only 

3 for  precip  probe  only 

4 for  VCO  data  only 
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ex.  to  dv  up  decimally,  2 files  every  10th  record  for 
Precj  o probe 

-DEC--2-10 3 - = blank 
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2.1.4  Output  Details 

The  DEC  (decimal)  and  SEL  (selected)  output  listings 
are  shown  in  figure  1A  and  are  identical  in  format.  The 
R=  and  F=  indicate  the  record  and  file  numbers  respectively. 
The  format  of  this  output  is  the  same  as  the  64  word  record 
structure  illustrated  in  Appendix  3. 

The  probe  select  option  of  the  decimal  output  is  shown 
in  figure  IB . ^ In  this  case  the  VCO  data  was  chosen  (a  4 in 
column  15  of  the  DEC  card) . The  first  number  on  each  data 
line  is  the  PMS  elapsed  second  counter  followed  by  the  probe 
designator,  in  this  case  VCO  data.  The  thirteen  VCO  values 
follow.  The  real  time  clock  (hours,  minutes  and  seconds) 
output  is  at  the  end  of  the  line. 

The  last  set  of  output  (figure  1C)  is  a listing  of  the 
status  parameters.  The  literal  BO,  CO  or  DO  denote  the  probe 
scatter,  cloud  or  precip.  The  first  column  lists  the  elapsed 
second  counter.  The  next  column  shows  the  real  time  clock. 
The  number  in  parenthesis  at  the  top  of  this  column  is  the 
elapsed  second  corresponding  to  this  time.  Thus  since  this 
number  is  2,  the  elapsed  seconds  corresponding  to  this  time 
is  two  more  than  that  given  in  the  first  column.  The  other 
columns  are  the  status  values  corresponding  to  the  elapsed 
second  counter  ending  in  that  particular  digit.  Refer  to 
Appendix  1 for  an  explanation  of  these  housekeeping  values 
and  to  Appendix  4 for  a listing  of  the  parameters  monitored 
each  second. 
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2.2.1  Program  Description 

Program  KNOLLlD  is  the  primary  post  processing  program 
for  data  collected  by  the  PMS-lD  systems.  These  devices  are 
currently  on-board  the  C130E  and  LEAR  aircraft.  At  an  earlier 

time,  however,  data  were  also  collected  from  a similar  device 
on  an  older  C130A  aircraft. 

Program  KNOLLlD  provides  the  user  with  a complete  parti- 
cle density  distribution  spectra  in  the  range  of  1 to  4650 
microns.  From  this  distribution  a calculated  set  of  meteoro- 
logical parameters  (i.e.  Liquid  Water  Content  (M) , Radar  Re- 
flectivity (Z) , Median  Volume  Diameter  (Dq) , Spectral  Parameter  (K) , 
Form  Factor  (F) , and  Stability  (S) , are  produced  in  tabular 
form.  Environmental  parameters  such  as  temperature,  dewpoint, 
etc.  are  also  calculated  and  listed.  In  addition  M,  Z and  D 

o 

are  produced  graphically  as  line  printer  plots. 

KNOLLlD  is  the  first  program  executed  in  the  data  proces- 
sing stream.  It  provides  processed  values  in  the  form  of  tapes 
and  cards  for  additional  production  programs. 

The  program  is  run  under  control  of  the  CDC  6600  batch 
processor.  Input  is  via  both  cards  and  tape.  The  Cl 30  tape 
is  recorded  on  a 7-track  Kennedy  recorder;  the  LEAR  system 
uses  a nine  track  Pertec  recorder  with  a slightly  different 
block  size.  Because  of  this  the  LEAR  tape  must  be  pre-processed 
using  program  LEARPMS  (see  Section  2.9).  The  output  consists 
of  listings,  line  printer  plots,  tapes  and/or  cards. 
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FILE  DESC.  IPTION 


UNIT 

TYPE 

MEDIUM 

1 

INPUT 

TAPE 

2 

OUTPUT 

TAPE 

3 

OUTPUT 

PAPER 

4 

OUTPUT 

CARD 

5 

INPUT 

CARD 

6 

OUTPUT 

PAPER 

7 

OUTPUT 

TAPE 

9 

OUTPUT 

PAPER 

10 

OUTPUT 

PAPER 

DESCRIPTION 

PMS-1D  TAPE 
OUTPUT  SUMMARY  TAPE 
LISTING  OF  M,  Z,  Dq 
M AND  Z 
INPUT  DATA 
ADDITIONAL  OUTPUT* 

M AND  Z (USED  IN  RAPP) 

LINE  PRINTER  PLOT  OF  M,  Z,  Dq 
OUTPUT  FILE 


* The  output  file  (TAPE10)  may  either  be  the  default  ABRES 
long  listing  or  the  abbreviated  AFWL  listing.  In  either  case 
TAPE6  will  contain  the  output  NOT  selected  and  may  be  viewed 
by  adding  the  appropriate  control  cards. 
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2.2.2  Control  Cards 

JOBNR , CM66000 , T600 , TP1 . PROB  NO.  NAME 

PAUSE.  PLS  MOUNT  DISK  LYCPFI 
MOUNT , SN=LYCPFI , VSN=LYCPF I . 

SETNAME (LYCPFI) 

ATTACH, LG O, KNOLLlDBIN , ID=GLASS , MR= 1 . 

REQUEST  TAPEl , S , HI , VSN=PMSXXX . 

REQUEST, TAPE2 , *PF , SN=LYCPFI . * 

REQUEST, TAPE7, *PF, SN= LYCPFI . * * 

FILE  (TAPEl, RT=U , BT=K,MRL=1024 , MBL=1024 , RB= 1 , BFS=10  5 ) 
LDSET , FILES=TAPEl , PRESET^  ZERO . 

MAP, OFF. 

LGO . 

CATALOG , TAPE 2 , PLOTDATA, ID=GLASS . * 

CATALOG, TAPE 7, RAPPDATA, ID=GLASS. ** 

REWIND , TAPE  3 , TAPE  8 , TAPE  9 . 

COPY, TAPE 3. 

COPY, TAPE 8. 

COPY, TAPE 9. 

REWIND, TAPE6. *** 

COPY , TAPE 6 . * * * 

EXIT (S ) 

REWIND, TAPE6. *** 

COPY , TAPE 6 . * * * 

CATALOG , TAP  E 2 , P LOTDATA , I D=  GL AS  S . * 

CATALOG, TAPE 7 , RAPPDATA, ID=GLASS . * * 

REWIND, TAPE3, TAPE8 ,TAPE9 . 

COPY, TAPE 3. 

COPY, TAPE 8. 

COPY , TAPE9 . 

7/8/9 


26 


DATA  CARD 

7/8/9 

6/7/S/9 


k 

k * 
k k k 


for  Plottape 
for  Rapptape 
for  optional 


option  must  have  REQUEST  TAPE2  card 
option  must  have  REQUEST  TAPE 7 card 
output 
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2.2.3  Data 

card 


card 

cards 

card 


Cards 

1:  ID  CARD 


cc  1-10 

FLT  E78 

-NN 

11-20 

DD  MON 

YY 

21-30 

PMSXXX 

(INPUT 

TAPE) 

31-40 

LYCXXX 

(PLOT 

TAPE ) 

41-50 

LYCXXX 

(RAPP 

TAPE) 

TIME  CARD 

cc  1-15 

START: 

HH : MM : SS 

22-35 

STOP: 

HH : MM: SS 

43-65  ZERO  SECONDS  = HH:MM:SS 


4:  NAMELIST 

’ CARDS 

$ SCOEF  $ 

(SOUNDING  COEFFICIENTS) 

$ VCOEF  $ 

(VCO  CALIBRATIONS) 

OPTION  CARD 

cc  1-5 

ICRYS 

6-10 

ICLK 

11-15 

I DAT 

20 

IPLT 

21-25 

INTA 

30 

ITMP 

31-35 

JVCO 

36-40 

NSKP 

45 

IDEK 

28 


50 

I RAD 

55 

INTRP 

60 

IDMZ 

65 

IVEL 

70 

IFORM 

card  6:  OPTIONAL  DATA 

any  RADIOSONDE  or  JVCO  data  required  by  switches 
set  on  the  option  card 

card  7 : 


The  following  cards  7A-7D  can  be  interspersed 


card  7A: 


card  7B: 


TYPE  CARDS 

(OPTIONAL  - 

15  MAXIMUM) 

1-4 

EDIT 

6-13 

HH : MM: SS 

(START) 

16-23 

HH : MM: SS 

(STOP) 

25-26 

CLOUD  TYPE 

27-28 

PRECIP  TYPE 

EDIT  CARDS 

(OPTIONAL  - 

5 MAXIMUM) 

1-4 

EDIT 

6-13 

HH : MM: SS 

(START) 

16-23 

HH : MM: SS 

(STOP) 

26 

PROBE  1 

28 

PROBE  2 

30 

PROBE  3 

31-54 

CHANNELS  TO 
8 MAXIMUM) 

BE  EDITED  (13  FORMAT 

29 


card  7C:  HTOX  CARDS  (OPTIONAL  - NO  MAXIMUM) 


1-4  HTOX 

6-7  PARTICLE  TYPE 

9-10  EQUATION  NUMBER 

12-13  ARGUMENT  TO  BE  CHANGED  (l=m,  2 = b, 

3=BREAKPOINT (N) ) 

15-30  NEW  VALUE  (REAL  - F16.0) 


card  7D : XTOD  CARDS  (OPTIONAL  - NO  MAXIMUM) 


1-4 

XTOD 

6-7 

PARTICLE  TYPE 

9-10 

EQUATION  NUMBER 

12-13 

ARGUMENT  TO  BE  CHANGED  (l=co 

3= BREAKPOINT (C) ) 

15-30 

NEW  VALUE  (REAL  - F16.0) 

30 


NAMELIST  DESCRIPTION 


card  3:  JSCOEF  card 


This  card  allows  the  coefficients  of  the  pressure- 
height  sounding  equation  to  be  changed. 

Control  variable:  S(I)  where  I = 1 to  5 

1^ 

1 first  order  coefficient 

2 second  order  coefficient 

3 third  order  coefficient 

4 fourth  order  coefficient 

5 fifth  order  coefficient 

For  example,  to  change  the  second  order  coefficient, 
the  control  variable  is  S(2)  = new  value. 


card  4 : $VCOEF  card 


This  card  allows  the  VCO  calibration  coefficients 
to  be  input.  The  default  coefficients  are  0,1  and 
0.  Thus  each  run  must  have  the  desired  coefficients 
input. 

Control  variable:  C(I,J) 

where  I = 1 to  3 and  J = 1,10  for  E Model 

J = 1,5  for  LEAR 


J E MOD 


LEAR 


1 intercept 

2 slope 

3 third  order  coef. 


1 

2 

3 

4 

5 

6 


indicated  airspeed 

temperature 

EWER 

n/a 

dew  t 
LWc  “u * 


air  pressure 
delta -pressure 
temperature 
dewpoint 
LWC- JW 
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J E MOD 

7 heading 

8 air  pressure  (Kistler) 

9 true  airspeed 
10-13  n/a 


For  example,  to  input  the  third  order  coefficient 
of  JW,  the  control  variable  is  C(3,6)  = value. 


card  5:  OPTION  CARD  DESCRIPTION 


cc  5 

10 

15 

20 

22-25 

30 

35 

40 


ICRYS  = 1 
= 2 
= 3 

ICLK  = 1 
= 2 

IDAT  = 0 
> 0 

IPLT  = 0 
= 1 

INTA  = 0 
> 0 

ITMP  = 0 
= 1 
= 2 

JVCO  = 0 
> 0 

NSKP  = 0 
> 0 

< 0 


TYPING  BY  TEMPERATURE 
UNUSED 

TYPING  BY  MANUAL  INPUT 
TIME  FROM  AIRCRAFT  CLOCK 
TIME  FROM  PMS  BUFFER 
PRELIMINARY  DATA 
FINAL  DATA 

PLOT  TAPE  NOT  PRODUCED 
PLOT  TAPE  PRODUCED 
FULL  DATA 

AVERAGING  INTERVAL  IN  SECONDS 

TEMPERATURE  DETERMINATION 
BY  VCO 

DETERMINATION  BY  STANDARD 
ATMOSPHERE 

TEMPERATURE  DETERMINATION  BY 
RADIOSONDE  PROFILE 

CALCULATE  VCO  PARAMETERS 

NUMBER  OF  JVCO  CARDS 

PROCESS  FIRST  FILE 

NUMBER  OF  FILES  TO  SKIP  BEFORE 
PROCESSING 

NUMBER  OF  RECORDS  TO  SKIP  BE- 
FORE PROCESSING 
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• C 4 5 

IDEK 

= 

0 

DECK  NOT  PUNCHED 

= 

1 

DECK  PUNCHED 

50 

I RAD 

= 

0 

RAPP  TAPE  NOT  PRODUCED 

= 

1 

RAPP  TAPE  PRODUCED 

55 

INTRP 

= 

0 

NO  DATA  INTERPOLATION 

= 

1 

DATA  INTERPOLATION 

60 

IDMZ 

= 

0 

NO  DMZ  PLOT 

= 

1 

DMZ  ■ PLOT 

65 

IVEL 

= 

0 

USE  TAS  VCO 

= 

1 

USE  STANDARD  CALCULATION 

70 

IFORM 

= 

0 

ABRES  OUTPUT  TO  FILE10 

= 

1 

AFWL  OUTPUT  TO  FILE10 

NOTES 

1.  ICRYS 

For  ICRYS=3  at  least  1 TYPE  card  is  required 


2 . ICLK 

For  ICLK=1  PMS  on  time  may  be  set  to  zero 

3.  IDAT 

for  I DAT > 0 IDAT  is  printed  as  FINAL  # on  page  2 of 
output 


4. 


ITMP 


For  ITMP= 1 standard  atmosphere  equation  is: 


T 


76.88288 (P) 


0.190284 


For  ITMP=2  profile  uses  a maximum  of  20  levels 


5 . JVCO 

For  JVCO>0  allows  VCO  parameters  to  be  manually  input 
for  specified  time.  This  option  overrides  the  ITMP 
selection 


JVCO  CARDS  DESCRIPTION  (OPTIONAL) 


JVCO  cards  are  required  only  if  cc  34-35  on  the  option 
card  is  greater  than  zero.  These  cards  override  the  VCO 
calculations;  they  supply  the  basic  meteorological  parameters 
necessary  for  processing.  The  input  values  are  constant  for 
the  duration  of  the  time  specified.  The  card  format  follows. 


cc 

FORMAT 

DESCRIPTION 

1-10 

A10 

@HH:MM:SS@  interval 

start  time 

11-20 

A10 

@HH : MM: SS@  interval 

stop  time 

21-30 

F10 . 0 

PRESSURE  (mb) 

i 

free  field 

31-40 

F10 . 0 

AIRSPEED  (knots)  1 

- 1 

with 

41-50 

F10.0 

TEMPERATURE  (°C) 

decimal 

51-60 

F10 . 0 

HEIGHT  (feet) 

point 

RADIOSONDE  CARDS  DESCRIPTION  (OPTIONAL) 

These  cards  are  required  only  if  the  option  card  cc  30 
is  2.  They  override  the  VCO  temperature  calculation  by  al- 
lowing the  temperature  to  be  determined  by  linear  interpola- 
tion using  the  RADIOSONDE  temperature-pressure  profile.  The 
first  card  contains  the  number  of  cards  to  follow  (20  maximum 
The  remaining  cards  contain  the  RADIOSONDE  profile.  These 
cards  must  be  in  descending  pressure  sequence. 


card  1 


cards  2 . 


cc 

1-2 

number  of  levels 

.n 

cc 

1-10 

PRESSURE  (mb) 

11-20 

TEMPERATURE  (°C) 
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TYPE  CARDS  DESCRIPTION  (OPTIONAL) 

The  ype  cards  are  required  only  if  the  option  card 
cc  5 is  3.  These  cards  specify  the  particle  type  to  be 
used  for  the  selected  interval.  The  cloud  probe  type  and 
precip  probe  type  may  be  different.  The  scatter  probe  is 
always  type  rain. 

TYPE  CARDS 


cc 

FORMAT 

DESCRIPTION 

1-4 

A4 

TYPE 

5-14 

A10 

START  TIME 

(@HR: MN : SC@  ) 

15-24 

A10 

STOP  TIME 

( @HR: MN : SC@ ) 

25-26 

12 

TYPE  KEY* 

(CLOUD  PROBE) 

27-28 

12 

TYPE  KEY* 

(PRECIP  PROBE) 

* The 

TYPE  KEY  used 

is  obtained  from 

the  particle  type 

table 

on  page 
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EDIT  CARDS  DESCRIPTION  (OPTIONAL) 


These  cards  allow  selected  particle  channel  counts  to 
be  edited.  The  current  editing  constraints  are  as  follows 


To  edit  channel  1 or  15 
To  edit  channels  2-14 
To  edit  channels  1 & 2 


- there  must  be  two  valid  channel 
on  one  side 

- there  must  be  one  valid  channel 
on  both  sides 

- there  must  be  a valid  channel  3 


cc 

FORMAT 

DESCRIPTION 

1-4 

A4 

EDIT 

5-14 

A10 

START  TIME 

( @HR: MN : SC@ ) 

15-24 

A10 

STOP  TIME 

(@HR:MN : SC@ ) 

25-26 

12 

PROBE  NO.  1 

i / 1 = scatter 

27-28 

12 

PROBE  NO.  1 

l J 2 - cloud 

29-30 

12 

PROBE  NO . J 

I 3 = precip 

31-33 

13 

CHANNEL  NO 

34-36 

13 

CHANNEL  NO 

37-39 

13 

CHANNEL  NO 

40-42 

13 

CHANNEL  NO 

43-45 

13 

CHANNEL  NO 

46-48 

13 

CHANNEL  NO 

49-51 

13 

CHANNEL  NO 

52-54 

13 

CHANNEL  NO 
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HTOX  CARDS  DESCRIPTION  (OPTIONAL) 

The  TOX  cards  are  used  to  change  the  default  intercept 
(b) , slope  (m)  or  breakpoint  (N)  parameters  for  each  adjusted 
channel  size  equation.  For  example,  to  change  the  second 
slope  for  particle  type  17  (Aggregate  Plates  and  Dendrites) 
from  1.20  to  1.25  the  following  HTOX  card  is  required: 


cc 

VALUE 

DESCRIPTION 

1-4 

HTOX 

CARD  CODE  (FIXED) 

6-7 

17 

PARTICLE  TYPE 

9-10 

02 

SECOND  EQUATION 

12-13 

01 

SLOPE 

15-30 

1.25 

NEW  VALUE 

XTOD  CARDS 

DESCRIPTION  (OPTIONAL) 

The  XTOD 

cards  are  used  to  change  the  default 

coefficient 

exponent  (ex) 

, or  breakpoint  (C)  for  each  equivalent  melted 

diameter  equation.  For  example,  to  change  the  third  exponent 

for  particle 

type  2 3 

(Rimed  Dendrites)  from  1.0  to 

0.9  the 

following  XTOD  card  is  required: 

cc 

VALUE 

DESCRIPTION 

1-4 

XTOD 

CARD  CODE  (FIXED) 

6-7 

23 

PARTICLE  TYPE 

9-10 

03 

THIRD  EQUATION 

12-13 

02 

EXPONENT 

15-30 

0.9 

NEW  VALUE 
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PARTICLE  TYPE  TABLE 


TYPE 


NAME 


1 

RAIN 

3 

WET  SNOW 

5 

LARGE  SNOW 

7 

SMALL  SNOW 

9 

BULLET- ROSETTES 

11 

COLUMNS 

13 

NEEDLES 

15 

PLATE  F 

17 

AGG  P & D 

19 

DENDRITE  FAMILY 

21 

GRAUPEL 

23 

RIMED  DENDRITE 

25 

NOT  USED 

27 

NOT  USED 

29 

NOT  USED 
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.2.4  Output  De  ails 

The  , rst  five  pages  of  a KNOLLlD  printout  show  the 
input  and  output  parameters  for  this  particular  run  (fig.  2A-2E) . 

Page  one  consists  of  a "banner"  page  followed  by  the 
version  number  of  the  KNOLLlD  program.  Page  two  (figure  2B) 
is  described  as  follows 

1)  flight  identification  and  flight  date 

2)  date  and  time  of  processing 

3)  the  time  corresponding  to  zero  elapsed  seconds  also  which 
clock  is  to  be  used  "A/C"  or  "PMS" 

4)  the  start  and  stop  times 

5)  the  input  tape  number  (blank  in  this  case) 

6)  if  positive  the  number  of  files  to  skip  before  processing, 
if  negative  the  number  of  records  (0  in  this  case) 

7)  averaging  interval  (4) 

8)  whether  interpolation  is  to  be  used  or  not 

9)  whether  a plot  tape  is  to  be  made  and  what  its  number  is 

10)  whether  a RAPP  tape  is  to  be  made  and  what  its  number  is 

11)  whether  the  line  printer  plot  of  M,  Z and  Dq  is  to  be  made 

12)  whether  punched  card  output  is  to  be  made 

13)  the  form  of  the  output,  preliminary  or  final 

14)  method  of  temperature  calculation;  either  by  VCO  or 
radiosonde  profile 

15)  how  particle  typing  is  determined;  by  temperature  or  type 
cards.  If  the  type  card  option  is  used;  a listing  of 
probe,  type  and  effective  times  is  shown. 

Page  three  lists  the  VCO  calibrations  and  probe  configura- 
tions used  for  this  particular  run.  If  the  values  for  these 
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VCO 1 s should  be  changed  see  the  $VCOEF$  card. 

Page  four  lists  the  adjusted  class  size  equations  for 
each  particle  type,  this  table  may  be  changed  by  using  IITOX 
cards . 

Page  five  lists  the  equivalent  melted  diamter  equations 
for  each  particle  type,  this  table  may  be  changec  by  using 
XTOD  cards. 

There  are  two  different  output  formats  generated  by 
KNOLLlD;  the  standard  ABRES  (Advanced  Ballistic  Re-Entry 
Systems)  and  the  abbreviated  AFWL  (Air  Force  Weapons  Lab.) . 
They  both  have  two  average  intervals  per  output  page.  One 
of  the  ABRES  intervals  will  be  explained  here.  This  will 
facilitate  the  AFWL  description.  All  references  are  to 
figure  2F. 

A)  The  number  of  one  second  data  samples  that  were  averaged 
to  make  this  table 

B)  Which  clock  was  used  to  derive  time  "AIRCRAFT"  or  "PM S" 

C)  the  start  and  stop  time  of  this  interval 

D)  The  time  as  calculated  by  converting  elapsed  seconds  (I) 
to  hours,  minutes,  and  seconds  and  adding  this  to  the  PMS 
on  time  (H) 

E)  Flight  identification 

F)  Particle  typing  indicator 

G)  Data  analysis  indicator 

H)  The  real  time  that  corresponds  with  a PMS  TIME  COUNT  = 0 

I)  The  number  of  seconds  counted  by  the  PMS-lD  hardware  since 
it  was  turned  on.  (three  digits  on  the  data  tape  - KNOLLlD 
adds  the  leftmost  fourth  digit  upon  rollover) 
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J)  The  cl  .nnol  number  for  reference 

K)  The  c(  ter  diameter  for  this  channel  of  this  probe.  It 
is  caj  'lated  by  adjusting  the  channel  number  by  HTOX 
equations  (per  particle  type)  and  then  being  melted  down 
by  XTOD  or  equivalent  melted  diameter  equations 

L.)  The  normalized  number  density.  It  is  calculated  by  com- 
puting the  number  of  particles  that  would  be  detected  by 
this  channel  size  in  a cubic  meter  of  sample  volume.  Then 
for  comparison  with  the  other  channels  it  is  normalized 
by  dividing  it  by  the  channel  barwidth. 

K and  L above  are  repeated  once  for  each  channel  of  each  probe. 

Scatter  on  the  left,  cloud  in  the  middle,  and  the  precip  probe 

on  the  right. 

M)  This  is  the  column  of  calibrated  VCO  values.  FRSOTPOINT 
becomes  DEWPOINT  when  the  temperature  is  greater  than  or 
equal  to  0°C.  There  are  two  type  of  airspeed,  CALC  as 
derived  from  pressure,  A pressure  and  temperature,  or  TRUE 
which  is  given  by  the  true  airspeed  - computer . The  asterisk 
next  to  (m/sec)  shows  which  airspeed  is  used  as  velocity 

in  the  calculation  of  sample  volume.  It  may  be  changed 
by  changing  the  IVEL  variable  on  the  option  card.  The 
values  of  "N/A"  for  a particular  VCO  means  it  is  not  appli- 
cable to  this  aircraft  or  version. 

N)  Total  particle  count  in  this  size  channel 

O)  The  calculated  liquid  water  content  (LWC)  of  each  probe 
is  listed  along  with  the  cumulative  average  (from  the 
beginning  of  this  pass  - as  determined  by  TYPE  cards)  in 
parenthesis 

P)  The  radar  reflectivity  (Z)  of  each  probe. 
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Q)  The  spectral  parameter  (K)  (i.e.  the  ratio  of  LWC  to  square 
root  of  Z)  for  each  probe 

R)  The  median  volume  diameter  (Dq)  for  each  probe. 

S)  If  two  TYPE  cards  ccntr ibuted  to  this  average  this  would 
read  "CHANGED"  rather  than  "FIXED"  and  if  that  were  the 
case  this  particular  average  would  be  useless  since  the 
the  values  were  calculated  with  different  diameters  and 
barwidths . 

T)  The  particle  types  used  to  calculate  center  diameters 

U)  FORM  FACTOR  a non  physical  parameter  FF  = K/C/N^  where 

C = g-  P10  , and  NT  is  the  sum  of  number  density  for 

channel  five  of  the  cloud  probe  up  to  channel  15  of  the 
precip  probe 

V)  Stability  factor  (s)  a number  which  reflects  the  homo- 
geneity of  the  sample  S = /N  /FF 

AA)  Totals  (0-R)  for  scatter  probe 

BB)  Totals  (0-R)  for  cloud  probe 

CC)  Totals  (0-R)  for  precip  probe 

DD)  Totals  (0-R)  for  the  cloud  and  precip 

EE)  All  channels  marked  with  a # are  channels  whose  sampling 
area  overlaps  that  of  the  next  probe.  In  this  case  chan- 
nels 11-15  of  the  cloud  probe  overlaps  channel  one  of  the 
precip,  so  these  channels  are  not  included  in  the  totals 
of  cloud  and  precip. 

The  AFWL  output  is  very  similar  to  the  ABRES  in  content, 
however  the  format  is  different.  It  is  designed  to  fit  within 
an  8 1/2  x 11  paper  and  thus  two  "!"  are  printed,  one  at  the 
top  of  the  page  and  one  at  the  bottom  to  form  a cutting  refer- 
ence line.  The  significant  differences  in  calculations  between 
the  AFWL  and  ABRES  output  are  explained  below,  see  figure  2G. 
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1)  Sample  <j  start  time 

2)  Only  1 ie  precip  probe  type  is  listed 

3)  The  ct  ter  diameter  size,  it  is  only  adjusted  and  not 
melted 

4)  The  normalized  number  density  is  calculated  by  dividing 
through  by  the  unmelted  barwidths 

5)  The  LWC  in  gm/M3  for  each  probe 

6)  The  (0Q)  median  volume  diameter  for  each  probe 

7)  P - pressure  in  millibars 

8)  ALT  - altitude  in  kilometers 

9)  TEMP  - true  temperature  in  centigrade 

10)  FROSTPOINT/DEWPOINT  - in  centrigrade 

11)  TAS  - the  selected  airspeed  in  meters  per  second 

12)  Since  LWC  and  Dq  are  both  calculated  from  the  melted 
equations,  these  totals  reflect  the  elimination  of  any 
overlap  between  the  cloud  and  precip  probes. 

File  TAPE 3 in  figure  21  is  self  explanatory. 

File  TAPE 9 (figure  2H)  is  the  line  printer  plot  of  total 
D(  , LWC  and  Z versus  time.  The  plot  is  on  log  scale  with  M 
indicating  LWC,  Z for  Z and  D representing  D 
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2.3  KNPLT1D 

2.3.1  Program  I .cription 

This  program  produces  six  different  types  of  plots  dis- 
playing the  processed  PMS-lD  data  in  different  forms.  The 
input  tape  used  is  produced  by  program  KNOLL ID  (see  Appendix 
5 for  the  tape  format);  it  contains  averaged  data,  one  record 
per  average  interval.  In  addition  this  program  can  read  a 
processed  tape  containing  LEARJET  data.  This  tape  is  processed 
by  another  contractor  and  sent  to  LYC.  The  format  for  this 
tape  is  shown  in  Appendix  10. 

Each  plot  type  to  be  produced  requires  one  plot  request 
card.  This  card  contains  the  information  necessary  to  produce 
the  plot  (start  time,  stop  time,  etc.)  with  the  appropriate 
title.  The  details  of  the  plot  request  cards  are  included  in 
Section  2.3.3. 

To  reduce  plotting  time  and  paper  consumption  the  plots 
are  produced  on  35  mm  film.  The  following  sections  explain 
in  detail  the  plots  available. 


Control  Cards;  (LEAR) 


5 3 


JOBNR, CM1 70000 , T600 , TPl , STMFB , NTl . PROB . NO . 

REQUEST, TAPE 3 9 , RING , MT, VSN=WTCRT . 

PAUSE.  PLS  PUT  WT  NUMBERS  IN  DAYFILE 
VSN , TAPEl=LEARTP/NT . 

PAUSE.  TAPE 1 IS  9 TRACK  PE 
REQUEST , TAP El , PE , L . 

PAUSE.  PLS  MOUNT  DISK  LYCPFI 
MOUNT , SN=LYCPFI , VSN=LYCPFI . 

ATTACH, CRT, OFFLINECRT. 

SETNAME (LYCPFI ) 

FILE ( TAPE 1 , RT=U , BT=K , MRL=2 100, MBL-2 100, RB= 1 , BFS=2 12 ) 
LIBRARY (CRT) 

ATTACH , KP , KNPLT1D, ID=GLASS , MR  1 . 

ATTACH , P0 , PLOTLIBO , I D= GLASS , MR=1 . 

ATTACH , P 1 , P LOTL I B 1 , I D=  GL AS  S , M R= 1 . 

ATTACH , P 3 , PLOTLIB3 , ID=GLASS , MR= ] . 

ATTACH, P4 , PLOTLIB4 , ID=GLASS, MR  1 . 

LDSET , PRES ET=  ZERO , FILES -TAPE  1 . 

LOAD ( KP ,P0,P1,P3,P4) 

EXECUTE. 

7/8/9 

-DATA  CARDS- 

6/7/8/9 


NAMJ 


2.3.2  Control  (ids  (C130E) 


r>4 


JOBNR , CM  1 ;000 , T600,  TP 2 , STMFI3 . PROB.  NO. 

REQUEST, TAPE 39 , U 1 NO , MT, VSN  WTCRT . 

PAUSE.  PUS  PUT  WT  NUMBERS  IN  DAYFILE 
PAUSE.  PUS  MOUNT  DISK  UYCPFI 
MOUNT  , SN=UYCPFI , VSN -UYCPFI  . 

ATTACH , CRT , OFFU I NECRT . 

SETNAME (UYCPFI ) 

UIBRARY (CRT) 

REQUEST, TAPE2,MT,VSN=LYCXXX  (FROM  KNOUU1D) 

ATTACH , KP , KNPUTl D, ID=GUASS , MR=1 . 

ATTACH ,P0, PUOTUI BO , ID=GUASS, MR=1 . 

ATTACH , P 1 , PLOTUIB 1 , I D=GUASS , MR= 1 . 

ATTACH , P 3 , PUOTU L B 3 , I D=GUASS , MR= 1 . 

ATTACH , P4 , PUOTUI B4 , I D=GUASS , MR= 1 . 

UDSET, PRESET  ZERO 
LOAD ( KP , P 0 , P 1 , P 3 , P 4 ) 

EXECUTE. 

7/8/9 

-DATA  CARDS- 
6/7/S/9 


NAME 


2.3.3  Data  Cards 

The  first  data  card  is  the  information  card;  it  appears 
only  once. 

Each  plot  type  requires  a plot  request  card.  These  re- 
quest cards  are  unlimited  and  have  the  same  format.  The  cards 
are  divided  into  16  fields.  Each  plot  type  requires  certain 
fields  to  be  used;  all  the  unused  fields  may  be  left  blank. 


card  1 . 

information  card 

VAR 

cc 

FORMAT 

FUNCTION 

PLT 

1-3 

A3 

plot  type:  PEN  or  CRT 

CLK 

5 

11 

clock:  1=A/C  2=PMS 

IOUT 

7 

11 

0 = summary  only 

1 = date  & summary 

FLID 

11-20 

A10 

flight  id:  FLT  XYR-NN 

21-30 

A10 

date:  DD  MON  YR 

OPT 

4 5 

11 

0 = standard  data 

1 = LEARJET  data 

INT 

49-50 

12 

averaging  interval 

(OPT  = 1 only) 

cards  2...n  plot  request  cards 


cc 

FORMAT 

FUNCTION 

1 

11 

field  1 

3 

11 

field  2 

5-10 

16 

field  3 

15-20 

16 

field  4 

22-25 

14 

field  5 

cc 

FORMAT 

FUNCTION 

27-30 

14 

field  6 

32-33 

14 

field  7 

37-40 

14 

field  8 

42-43 

14 

field  9 

47-50 

14 

field  10 

52-55 

14 

field  11 

(Ji 

1 

O'* 

o 

14 

field  12 

62-65 

14 

field  13 

67-70 

14 

field  14 

72-75 

14 

field  15 

77-80 

14 

field  16 

Cards  2 through  n may  appear  in  any  sequence,  however,  con- 
siderable time  is  saved  if  cards  with  same  time  limits  are 
consecutive . 

The  following  page,  "Request  card  summary"  shows  the  required 
fields  for  each  plot  type. 


REQUEST  CARD  SUMMARY 


* 

* 


* 
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AXIS  codes  are  shown  in  Append! 
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SCATTER  PLOTS 

This  ption  produces  two  plots  of  LOG  M vs.  LOG  Z.  Plot 
one  uses  Precip  Probe  data  only  and  plot  two  uses  total  data. 
A least  squares  linear  fit  is  drawn  through  the  data.  The 
coefficients  appear  in  the  form  M = A*Z**B. 


1 

start  time  (HHMMSS) 
stop  time  (HHMMSS) 
pass  (rj) 

11TKM*  10  (rj) 

i.e.  30  km  = 300 
27.2km  = 272 

HISTOGRAMS 

This  option  produces  two  histograms,  the  first,  LOG  M 
vs.  time,  and  the  second,  LOG  Z vs.  time.  Each  histogram 
contains  three  plots,  one  for  each  probe.  The  time  axis 
is  set  for  a 300  second  maximum,  however,  less  data  may  be 
plotted. 

INPUT  CARD: 

field  1 2 

field  3 start  time  (HHMMSS) 

field  4 duration  seconds  (rj) 

field  5 pass  (rj) 

HTKM*10  (rj) 


INPUT  CARD: 

field  1 
field  3 
field  4 
field  5 
field  6 


field  6 
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VCO ' s 


This  option  produces  three  plots  witli  two  VCO's  per  plot. 
The  six  VCO's  plotted  are  not  variable.  They  are:  Magnetic 
Heading,  LWC(JW),  Temperature,  Dewpoint,  Pressure,  and  Accel- 
eration. The  scales  for  Acceleration  and  1,WC(JW)  are  fixed 
at  + lg  and  -.1  to  +.8  gm/m3  respectively.  The  remaining 
VCO's  have  fixed  ranges  but  the  scale  slides  to  plot  as  many 
points  as  possible.  The  fixed  ranges  are:  Pressure  (30  mb) , 
Temperature  (12°C),  Dewpoint  (12°C) , and  Heading  (deg).  The 
time  axis  is  set  for  a 300  second  maximum,  however,  less  data 
may  be  plotted. 

INPUT  CARD: 

field  1 
field  3 
field  4 


This  option  produces  six  plots  of  number  density  vs. 
equivalent  melted  diameter.  Each  plot  uses  only  data  from 
the  specified  probe  for  six  liquid  water  content  bands.  The 
LWC  bands  are  variable  but  the  probe  must  be  the  same  for 
all  six  plots.  If  the  particle  type  remains  constant  for 
the  entire  interval,  an  average  line  is  drawn  through  the 
data.  Also,  the  cumulative  percent  mass  is  superimposed  on 
each  plot. 


3 

start  time  (HHMMSS) 
duration  seconds  (rj) 

DENSITY  SPECTRA 


INPUT  CAI  > 


f ieJ 

1 

4 

field 

2 

probe 

(l=sc,  2 

field 

3 

start  time 

( HHMMSS ) 

field 

4 

stop  time 

(HHMMSS) 

field 

5 

mIn  1 BL 

r j 

field 

6 

max 

rj 

field 

7 

mln  } BR 

r j 

field 

8 

max 

r j 

field 

9 

min  } ML 

r j 

field 

10 

max 

r j 

field 

11 

min  } MR 
max 

r j 

field 

12 

r j 

field 

13 

min  j TL 

rj 

field 

14 

max 

rj 

field 

15 

min 

r j 

field 

16 

) TR 

max 

r j 

The  min  and  max  values  are  the  lower  and  upper  limits 
of  each  LWC  hand.  The  LWC  limits  are  in  units  of  mg/M3  with 
the  10“ 3 exponent  omitted.  Since  there  are  four  columns  per 
limit,  the  absolute  range  of  limits  is  from  1 = .001  mg/M3  = 
1 mg/M3  to  9999  ~ 10.0  gm/M ’ = 10 4 mg/M3. 

The  two  letter  code  BL,  MR,  etc.  indicate  which  plot  on 
the  page  is  used.  The  plots  are  oriented  as  follows: 

TL  TR 

ML  MR 

BL  BR 


f)1 


Less  than  six  plots  may  be  utilized  by  leaving  the  appro- 
priate columns  on  the  input  card  blank. 

MED I AN  VOLUME  DIAMETER 

The  Median  Volume  Diameter  module  produces  seven  plots 
of  the  following  form. 


D 

o 

D 

o 

vs . 

vs . 

(Z/M) 1/3 
(Z/M) 1/3 

Precip  only 

Cloud  & Precip 

D 

o 

vs . 

K 

Precip  only 

D 

o 

vs . 

K 

Cloud  & Precip 

ND 

o 

4/m 

vs . 

D/D 

o 

Precip  only 

0 

a 

■<; 

Vm 

vs . 

D/Dq 

Cloud  only 

ND 

o 

Vm 

vs . 

D/Dq 

Cloud  & Precip 

INPUT: 


field 

1 

5 

field 

3 

start  time 

(HUMMUS) 

field 

4 

stop  time 

(HHMMSS) 

field 

5 

pass 

(r  j ) 

filed 

6 

I1TKM*  10 

(r  j) 

Thii 


VCO  PLOT 

iodule  produces  a maximum  of  four  plots  per  input 
card.  The  plots  are  any  VCO,  LWC,  or  X versus  any  other 
VCO,  LWC  or  X.  There  arc  two  types  of  plotting:  scatter  or 
line.  This  option  does  not  have  a 300  point  maximum,  hence 
the  plots  can  be  run  for  an  entire  flight. 

In  the  special  case  of  plotting  JW  Liquid  Water  Content 
vs.  Height  an  additional  option  is  provided.  This  allows  for 
all  the  data  points  to  be  shifted  by  some  reference  equation. 
The  information  for  this  adjustment  is  on  a -$ADJUST  card. 
This  card  must  immediately  follow  the  input  card  which  speci- 
fies the  JW-LWC  vs.  HEIGHT  plot. 

INPUT  CARD: 


field 

1 

6 

field 

2 

# number  of  plots  this  card 

(4  max) 

field 

3 

start  time  (HHMMSS) 

field 

4 

stop  time  (HHMMSS) 

f ield 

5 

horizontal  axis 

(rj) 

field 

6 

vertical  axis 

Plot 

1 

(rj) 

field 

7 

type 

(rj) 

field 

8 

horizontal  axis 

(rj) 

field 

9 

vertical  axis 

Plot 

2 

(rj) 

field 

10 

type 

(rj) 

f ield 

11 

horizontal  axis 

(r  j ) 

field 

12 

vertical  axis 

Plot 

3 

(r  j ) 

f ield 

13 

type 

(r  j ) 

field 

14 

horizontal  axis 

(r  j) 

f ield 

15 

vertical  axis 

Plot 

4 

(rj) 

field 

16 

type 

(r  j ) 

6 3 


JW-LWC  ADJUSTMENT  CARD 

If  a VCO  plot  request  card  specifics  a LWC-JW  vs.  HEIGHT 
plot,  (i.e.  h axis  code  =5,  v axis  code  = 2,  and  type  code 
= 0)  the  next  card  must  contain  the  adjustment  parameters. 

The  required  parameters  arc: 

L = number  of  levels  ( LO  maximum) 

XJ  = origin  of  the  level 
SL  = slope  of  the  level 

HT  = height  (meters)  at  the  top  of  the 
level 

The  card  uses  a standard  namelist  format  with  the  control 
variable  being  $ADJUST.  If  the  option  is  not  desired  the 
card  must  be 


$ADJUST  L=0 , $END 
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Figure  4A:  KNPLT1D  Sample  Plots 
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Figure  4D:  KNPLTlD  Semple  Plots 
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Figure  4E:  KNPLTlD  Sample  Plots 
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Figure  4F:  KNPLTlD  Sample  Plots 
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Figure  4H:  KNPLTlD  Sample  Plots 
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Figure  41:  KNPLT1D  Sample  Plots 
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Figure  4J:  KNPLTlD  Sample  Plots 
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2 . 4 DENPLOT 

2.4.1  Program  Description 

DENPLOT  produces  graphs  of  the  log  of  normalized  number 
density  vs.  particle  diameter  measured  during  a given  interval 
of  time  from  a user  specified  start  time.  It  is  run  from  the 
Tektronix  graphics  terminal  interactively. 

The  program  uses  a standard  KNOLL ID  (TAPE 2 ) plot  tape 
which  has  been  previously  cataloged  onto  the  disk  as  its 
input  tape. 
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2.4.2  Operating  Instructions 


To  use  DENPLOT  from  the  Tektronix  terminal. 


should  execute  the  following  commands: 


the  user 


LOGIN  , NAME  , PASSWORD , TTYN04- 
M,  PLS  MOUNT  DISK  LYCPFII 
MOUNT , SN=LYCPFI , VSN=LYCPFIi 
SETNAME (LYCPFI ) + 

ATTACH , TEK , TEKLIB , SN=SHAREDi 
LIBRARY (TEK) I 

ATTACH , LGO , DENPLOTBIN , ID=FRANCISI 
ATTACH , TAPEl , KNOLL 1DPLOTTAPE , ID=NAMEI 
CONNECT, INPUT! 

CONNECT, OUTPUT* 

LGOi 

[INPUT  RESPONSES] 
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2.4.3  Input  Responses 

After  the  LGO  command  has  been  executed,  the  program 
will  begin  asking  for  various  parameters: 

(1)  The  flight  ID 

(2)  The  date  of  the  flight 

(3)  Which  probe  to  begin  the  plot  with 

(4)  Which  probe  to  end  the  plot  with 

(5)  The  desired  starting  time 

(6)  Which  clock  to  take  the  time  from 

(7)  The  number  of  seconds  worth  of  data  to 
be  plotted 

For  example,  the  session  in  Figure  5 produces  a plot  of 
the  normalized  number  density  of  all  three  probes.  It  uses 
600  seconds  worth  of  data  starting  at  19:10:10.  It  uses  the 
PMS  clock  to  determine  the  time.  The  flight  took  place  on 
26  February  1978  and  its  flight  ID  was  E78-03.  The  graph 
produced  by  this  session  is  shown  in  Figure  6. 
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2.5  DENTIME 

2.5.1  Program  Description 

DENTIME  was  written  to  plot  the  log-^Q  of  number  density 
vs.  time.  It  is  usually  executed  interactively  from  either  a 
standard  Intercom  terminal  or  from  the  Tektronix  plotting 
terminal . 

The  program  uses  a standard  KN0LL1D  plot  tape  which  has 
been  previously  cataloged  onto  the  disk  as  its  input  tape. 
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2.5.2  Control  Commands 


In  order  to  run  DENTIME  from  the  Tektronix  plotting  ter- 
minals the  user  must  execute  the  following  commands: 


LOGIN , NAME , PASSWORD , TTYNO 
M, PLS  MOUNT  DISK  LYCPFI 
MOUNT, SN=LYCPFI , VSN=LYCPFI 
SETNAME (LYCPFI) 

ATTACH , LGO , DENTIMEBIN , I D= LEACH 
ATTACH , TAPE1 , KN0LL1DPL0TTAPE , ID=NAME 
ATTACH , LIB , TEKLIB , SN=SHARED 
LIBRARY (LIB) 

CONNECT, INPUT 
CONNECT, OUTPUT 
LGO 


In  order  to  run  DENTIME  from  a standard  Intercom  terminal, 
the  user  must  execute  the  following  commands: 

LOGIN , NAME , PASSWORD , TTYNO 
M, PLS  MOUNT  DISK  LYCPFI 
MOUNT , SN=LYCPFI , VSN=LYCPFI 
SETNAME (LYCPFI) 

ATTACH , LGO , DENTIMEBIN , ID=LEACI1 
ATTACH , TAPE1 , KNOLL 1DPL0TTAPE , ID=NAME 
ATTACH , LIB , SMALLONLINE , SN=SHARED 
LIBRARY (LIB) 

CONNECT, INPUT 
CONNECT, OUTPUT 
REQUEST , SMPLOT , *Q 
LGO 

[INPUT  DATA]  see  section  2.5.3 
DISPOSE, SMP LOT, SP 
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2.5.3  Input  Com.  ands 

After  the  LGO  command  is  executed  (at  either 
the  program  will  begin  asking  the  user  for  various 
parameters  concerning  the  type  of  plot  to  be  made. 

It  asks  for  the  following  parameters: 

(1)  type  of  terminal  being  used 

(2)  number  of  minutes  to  be  plotted 

(3)  which  probe  type  to  use 

(4)  type  of  density  to  be  plotted 

(5)  date  of  the  flight 

(6)  maximum  y-value  desired 

(7)  which  clock  time 

(8)  starting  time 


terminal ) 
input 
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2.5.4  Output  Details 


ENTER  TYPE  OF  TERMINAL 

(STANDARD  INTERCOM*  1*  TEKTRON  I X”  2)  _1_ 

ENTER  THE  NUMBER  OF  MINUTES 
TO  BE  PLOTTED  100 

ENTER  PROBE  TYPE 

(SCATTER*  \»  CLOUD*  2,  PRECI  P=  3/  TOTAL*  A)  __4_ 

ENTER  TYPE  OF  DENSITY  TO  BE  PLOTTED 
(UNNORMALIZED*  U NORMAL  I ZED*  2)  1 

ENTER  THE  DATE  OF  THE  FLIGHT  2fc  FEE  78 

ENTER  MAXIMUM  Y DESIRED  10 

ENTER  CLOCK  TO  USE  ( A/C*  \,  FM  S=  2)  2 

ENTER  START  TIME  HH  MM  SS  18  18  10 


END  OF  PLOTS. 

STOP 

. 634  CP  SECONDS  EXECUTION  TIME 
..  DISPOSED  SMPLOT,  SP 


Example  of  a run  of  DENTIME 
(user  responses  are  underlined) 

Figure  7:  DENTIME  Sample  Session 

The  session  recorded  in  the  following  sample  will  create 
a plot  of  the  log^Q  of  unnormalized  total  number  density  vs. 
time  for  100  minutes  starting  at  18:18:10.  The  data  is  being 
entered  on  a standard  Intercom  terminal,  the  time  is  taken 
from  the  PMS  clock,  and  the  flight  took  place  on  the  26th  of 
February,  1978. 
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2 . 6 VHPLOT 

2.6.1  Program  Description 

Program  VHPLOT  produces  calibrated  pen  plots  of  reflected 
PMS-lD  VCO  and  status  word  data  as  read  from  the  Kennedy  tape. 
Since  any  set  of  five  different  VCO/STATUS  values  can  be  plotted 
on  one  frame  and  one  frame  may  have  up  to  one  hour's  worth  of 
data,  this  program  may  be  used  to  visually  demonstrate  the 
error  of  one  piece  of  equipment  in  relation  to  others. 

Input  to  VHPLOT  is  via  cards  through  the  CDC  6600  batch 
processor  and  input  tape  is  the  same  as  the  one  used  in 

KNOLLlD  and  KNlUTIL  (Appendices  1-3) . Output  consists  of  both 
tabulated  data  and  CALCOMP  pen  plots.  Plot  output  may  also 
be  routed  to  the  Tektronix  graphics  terminal. 
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2.6.2  Control  C rds 

JOBNM , CM6  0000, T2001  ,TPl.  ACT.  NAME 

REQUEST , TAPEl , S , HI , VSN=PMSXXX . 

PAUSE.  PLS  MOUNT  DISK  LYCPFI 
MOUNT , SN^LYCPFI , VSN=LYCPFI . 

ATTACH ,LGO, VHPLOTBIN , ID=GLASS , MR=1 , SN=LYCPFI . 

ATTACH , PEN , ONLINEPEN , MR=1 . 

DISPOSE , PLOT, *OL . 

MAP, PART. 

FILE (TAPEl , RT=U , BT=K,MRL=1024 ,MBL=1024 , RB=1, BFS=105) 
LIBRARY (PEN) 

LDSET , FILES=TAPEl , PRESET= ZERO . 

LGO. 

7/8/9 

DATA  CARDS 
6/7/ 8/9 


x 


ALLOW  APPROXIMATELY  60  sec/hour/f rame 
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1-6 

8-13 

25 

29-30 

35 

39-40 

45 

49-50 

55 

59-60 

65 

69-70 


DESCRIPTION 

aircraft  model  A or  E 

$ CHANGES  CALIBRATION  CARDS* 

$ END  INSERTED  HERE  AS  REQUIRED 

for  n frames 

(MAXIMUM  ONE  HOUR  DATA  PER  FRAME) 
(CARDS  MAY  BE  IN  ANY  TIME  SEQUENCE) 

START  TIME  HHMMSS 

STOP  TIME  HHMMSS 

PLOT  ONE  TYPE** 

PLOT  ONE  CODE 

PLOT  TWO  TYPE** 

PLOT  TWO  CODE 

PLOT  THREE  TYPE** 

PLOT  THREE  CODE 

PLOT  FOUR  TYPE** 

PLOT  FOUR  CODE 

PLOT  FIVE  TYPE** 

PLOT  FIVE  CODE 


CALIBRATION  CONSTANTS  AND  AXIS  LIMITS  (IN  COUNTS)  CAN 
BE  CHANGED 

TYPE  SPECIFIES  WHETHER  A PLOT  WILL  BE  V (FOR  VCO)  OR 
H (FOR  HOUSEKEEPING) , THE  CODES  ARE  ON  THE  FOLLOWING 
3 PAGES 
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CALIBRATION  CARDS 


There  are  three  control  variables  that  may  be  used  with 
the  $CHANGES  namelist  input.  The  variables  shown  below  allow 
the  calibration  coefficients  and  axis  limits  to  be  changed  as 
needed. 


VCOA(I,J)  controls  any  changes  pertinent  to  the 
C130A  VCO's 

VCOE(I,J)  controls  any  changes  pertinent  to  the 
C130E  VCO's 

HSVL ( I , J ) controls  any  changes  pertinent  to  the 
housekeeping  data  for  either  aircraft 

1=1  specifies  calibration  intercept 
= 2 specifies  calibration  slope 
= 3 specifies  minimum  counts 
= 4 specifies  maximum  counts 


J = 


for  A model  VCO’s 
for  E model  VCO's 

for  housekeeping  data  - from  the  PMS-1D  status 
word  code 


1-13 

1-9 

1-30 
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PMS  ID  VCO  CODE 


C130A 

J 

C130E 

Pressure 

1 

A P 

A P 

2 

Temp 

Mag  Head 

3 

EWER 

Temp 

4 

UNUSED 

Event/Cloud 

5 

Dewp/1011 

LWC/JW 

6 

LWC/JW 

Rain 

7 

Mag  Head 

Tacan  Bearing 

8 

Pressure 

Tacan  Distance 

9 

true  airspeed 

Acceleration 

10 

Dewp/1011 

11 

Ice 

12 

Pitch 

13 

PMS  ID  STATUS  WORD  CODE 


msec  = PMS  elapsed  second  clock  modulo  10 
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2.6.4  Output  Description 

The  output  from  VHPLOT  consists  of  two  sections  (figures 
9 A & B) . Figure  9A  is  the  input  option  section.  The  first 
line  consists  of  this  particular  start  and  stop  time  as  well 
as  the  aircraft  that  produced  the  input  data.  The  next  part 
of  this  section  describes  the  particular  plots  selected.  In 
this  example  the  lowest  frame  to  be  plotted  (plot  1)  is  tem- 
perature (TEMP) , and  the  limits  of  calibration  are  from  0 to 
10,000  counts.  The  calibration  for  this  plot  is 

TEMP  = 0.0101 (COUNTS)  - 50.8054 

Figure  9B  shows  a calibrated  listing  of  each  plotted 
value.  The  first  two  columns  show  the  elapsed  seconds.  The 
columns  of  VCO  values  are  self  explanatory.  An  important 
note  is  that  since  any  particular  housekeeping  value  occurs 
once  every  ten  second  every  other  housekeeping  value  is  listed. 
However  the  one  selected  for  plotting  is  indicated  by  the  word 
"VALUE"  next  to  it. 

The  plot  output  (figure  10)  shows  how  each  of  the  plots 
are  offset  so  that  each  axis  limit  can  be  easily  read.  Any 
fluctuation  in  any  of  the  housekeeping  curves  is  a sure  in- 
dication of  some  sort  of  device  problem. 
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2 . 7 LWCTIME 

2.7.1  Program  Description 

LWCTIME  was  written  to  plot  the  log  of  liquid  water  con- 
tent vs.  time  (using  either  the  aircraft  or  PMS  clock)  for  a 
given  number  of  minutes  from  an  input  start  time.  It  can  pre- 
pare plots  from  either  a standard  INTERCOM  terminal  (where 
the  output  is  disposed  to  a plot  queue)  or  the  Tektronix  gra- 
phics terminal  (where  the  output  is  displayed  in  real  time) . 

LWCTIME  uses  the  standard  KNOLLlD  plot  tape  (see  Appendix 
5)  as  its  input  which  has  been  previously  cataloged  onto  the 
disk. 


lor. 


.7.2  Operatim  Instructions 


To  use  LWCTIME  from  a standard  INTERCOM  terminal,  the 
user  must  execute  the  following  commands: 

LOG IN, NAME, PASSWORD, TTYNO! 

M,  PLS  MOUNT  DISK  LYCPFI! 

MOUNT,  SN=LYCPFI,VSN=LYCPFI! 

SETNAME (LYCPFI ) 1 

ATTACH , LIB , SMALLONLINE , SN=SHARED! 

LIBRARY (LIB) ! 

REQUEST , SMPLOT , *Q! 

ATTACH , TAPE1 , KNOLL 1DPL0TTAPE , ID=NAME! 

ATTACH ,LGO , LWCTIMEBIN , ID=LEACH! 

LGO! 

[INPUT  RESPONSES] 

DISPOSE , SMPLOT , SP! 


To  use  LWCTIME  from  the  Tektronix  plotting  terminal,  the 
user  must  execute  the  following  commands: 

LOGIN , NAME , PASSOWRD , TTYNO! 

M,  PLS  MOUNT  DISK  LYCPFI! 

MOUNT, SN=LYCPFI ,VSN=LYCPFI! 

SETNAME (LYCPFI) ! 

ATTACH , LIB , TEKLIB , SN=SHAREDi 
LIBRARY (LIB) ! 

ATTACH, TAPE 1, KNOLL 1DPLOTTAPE , I D=NAME! 

ATTACH , LGO , LWCTIMEBIN , ID= LEACH! 

LGO! 

[INPUT  RESPONSES] 
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2.7.3  Input  Responses 

Regardless  of  which  type  of  terminal  is  being  use,  the 
program  will  begin  asking  questions  about  the  graph  to  be 
plotted  after  the  LGO  command  has  been  executed.  The  user 
must  input  the  following  parameters: 

(1)  Whether  a standard  INTERCOM  terminal  or 
the  Tektronix  plotting  terminal  is  being 
used 

(2)  The  number  of  minutes  to  be  plotted 

(3)  Which  probe's  liquid  water  content  (or 
their  total)  is  to  be  graphed 

(4)  The  date  of  the  flight 

(5)  Whether  the  time  is  to  be  taken  from  the 
aircraft  clock  or  the  PMS  clock 

(6)  What  the  desired  starting  time  is 

For  example,  the  terminal  session  in-  Figure  11  will  pro- 
duce a graph  of  the  log  of  the  total  liquid  water  content  vs. 
time  for  100  minutes.  The  graph  will  start  at  18:18:10  and 
the  time  will  be  taken  from  the  PMS  clock.  The  date  of  the 
flight  was  February  26,  1978  and  the  plot  was  being  prepared 
from  a standard  INTERCOM  terminal.  Figure  12  shows  the  graph 
produced  as  a result  of  this  session. 
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Figure  11:  LWCTIME  Sample  Session 


LGO 

ENTER  TYPE  OF  TERMINAL 

< STAN DARD  IN  TERCOM3 \,  TEK TRON I X3  2>  1 

ENTER  THE  NUMBER  OF  MINUTES 
TO  BE  PLOTTED  100 

ENTER  PROBE  TYPE 

< SCATTER=  l*  CLOUD=  PRECI  P3  3..  TOTAL3  A)  A 
ENTER  THE  DATE  OF  THE  FLIGHT  26  FEB  78 
ENTER  CLOCK  TO  USE  C A/C=  FMS3  2)  2 
ENTER  START  TIME  HH  MM  SS  18  18  10 


END  OF  FLOTS. 

STOP 

1.196  CF  SECONDS  EXECUTION  TIME 
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2 . 8 DIAMETER 

2.8.1  Program  Description 

Program  DIAMETER  calculates  and  lists  the  equivalent 
center  diameters  and  channel  limits  for  each  particle  type 
on  either  the  C130E  or  LEARJET. 

DIAMETER  can  be  run  either  from  an  INTERCOM  terminal 
or  from  the  CDC  6600  batch  processor  with  cards.  The  only 
drawback  with  using  an  INTERCOM  terminal  is  the  lengthy  out- 
put, however  this  can  be  disposed  to  a high  speed  printer. 

There  is  no  input  required  and  the  only  output  is  a 
printout  which  has  one  particle  type  per  page.  There  are 
nine  particle  types  and  two  aircraft  resulting  in  eighteen 
pages  per  execution. 


Ill 


2.8.2  Control  cards 

JOBNM, CM65000 , T10 . ACT 

ATTACH , LGO , DIAMETERBIN , ID=FRANCIS . 

LGO. 

6/7/S/9 


INTERCOM  USAGE 


ATTACH, LGO, DIAMETERBIN, ID=FRANCIS. 
LGO. 

DISPOSE , OUTPUT , PR , IAC 


r 

* 

1 


NAME 


LCM  TV*C  1 MIN 

MMlNW  CHANNEL  IS  S HININUN  CHANNEL  IS 
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Figure  13A:  DIAMETER  Sample  Output 
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Figure  13B:  DIAMETER  Sample  Output 
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2 . 9 LEARPMS 

2.9.1  Program  Description 

Program  LEARPMS  reformats  the  PMS-1D  data  tape  produced 
by  the  Pertec  9 track  recorder  on  board  the  LEARJET.  The 
output  tape  appears  identical  to  that  produced  by  the  C130-E 
Kennedy  7 track  recorder.  LEARPMS  must  be  run  on  the  CDC  6600 
through  the  batch  processor  stream.  No  data  cards  are  needed. 

Input  and  output  consists  entirely  of  tapes.  The  input 
tape  formats  are  shown  in  Appendices  11  and  12.  The  output 
tape  format  is  shown  in  Appendix  2 . 
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2.9.2  Control  Cards 

JOBNM, 60000, T300,TP1,NT1.  ACC#  NAME 

PUASE . TAPE 9 IS  NINE  TRACK  CODED 
VSN=TAPE9=TAPENO/NT. 

VSN, TAPEl=TAPENO . 

REQUEST, TAPE9,HD,L, NT.  (9  TRACK  - NO  RING) 

REQUEST , TAPEl , S , HI , RING . 

ATTACH  LGO , LEARPMSBIN , ID=GLASS , MR=1 . 

FILE (TAPE 9 , RT=U , BT=K , MRL=5576 , MBL=5576 , RB=1 , BFS=560 ) 
FILE (TAPEl, RT=U, BT=K, MRL=1024 , MBL=1024 , RB=1 , BFS=105) 
LDSET , FILES=TAPEl/TAPE9 , PRESET= ZERO . 

LGO. 

7/8/9  (MULTI-PUNCH) 

6/7/8/9  (MULTI-PUNCH) 


2.9.3  Data  Cards 

none  required 

I 
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2.10  HIAC1D 

2.10.1  Program  Description 

Program  HIAC1D  summarizes  and  tabulates  the  Aeromet  pro- 
cessed data  generated  by  the  PMS-lD  Knollenberg  device  aboard 
the  LEARJET. 


HIAC1D  produces  two  output  tapes  (TAPE2  and  TAPE7)  similar 
to  those  written  by  KNOLL ID.  It  also  prints  a listing  which 
has  the  same  meteorological  parameters  as  KNOLLlD. 


The  input  to  HIAClD  is  a 9 track  data  tape  and  as  such, 
it  must  be  executed  from  cards.  Output  consists  of  a formatted 
listing,  a line  printer  plot,  a data  summary  and  two  output 
tapes.  The  following  is  a summary  of  I/O  units  used  in  HIAClD. 


UNIT 

1 

2 

3 

5 

7 

9 


DESCRIPTION 

input 

output  plot  tape 
output  summary 
card  input 
RAPP  tape  output 
line  printer  output 
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2.10.2  Control  Cards 

JOBNM, CM65000 , T200 ,NTl , TPO . * * ACT  # NAME 

PAUSE.  PLS  MOUNT  DISK  LYCPFI 
MOUNT,  SN=LYCPFI,VSN=LYCPFI. 

ATTACH , LGO , HIAC1DBIN , ID=FRANCIS , SN=LYCPFI . 

VSN , TAPE 1=LTAXXX/NT . 

REQUEST, TAPE 1, PE, L, NR, NT.  (9  TRACK  - NO  RING) 

REQUEST, TAPE 7, RING, MT, VSN=TAPENO. * (7  TRACK  - RING) 

REQUEST, TAPE2, RING, MT, VSN=TAPENO. * (7  TRACK  - RING) 

FILE (TAPEl , RT=U, BT=K,MRL=2100 , MBL=2 100 , RB=1 , BFS=211) 
LDSET , FILES=TAPEl, PRESET= ZERO . 

LGO. 

REWIND , TAPE 3 , TAPE 9 . 

COPY, TAPE 3. 

COPY, TAPE 9. 

EXIT (S ) 

REWIND, TAPE 3, TAPE9 . 

COPY , TAPE3 . 

COPY , TAPE9 . 

7/8/9 

DATA  CARDS 
6/7/8/9 


FOR  NO  RAPP  TAPE  REMOVE  REQUEST , TAPE7 ,.. . 

FOR  NO  PLOT  TAPE  REMOVE  REQUEST , TAPE2 ,.. . 

FOR  EACH  OUTPUT  TAPE  INCREMENT  TPO  to  TP1,  TO  TP 2 
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2.10.3  Data  Card' 


card  1 

col  1-6 
9-10 

15 


TAPE  INPUT  TAPE  NUMBER  (LTAXXX) 

NEOF  NUMBER  OF  END  OF  FILES  TO  SKIP 

BEFORE  PROCESSING 

PREREJ  =0  (DEFAULT)  USE  CHANNEL  1 PRECIP 

=1  INTERPOLATE  CHANNEL  1 PRECIP 


CARDS  2-N+l  (N  IS  NUMBER  OF  "PASSES") 


col 

2-9 

IBEG 

START  PASS  TIME 

HH:MM: SS 

12-19 

IEND 

STOP  PASS  TIME 

HH:MM: SS 

21-25 

I PASS 

PASS  NUMBER 

26-30 

INTA 

AVERAGING  INTERVAL 

33-65 

NED 

EDIT  PARAMETERS  (SEE  BELOW) 

NED 

EDIT 

PARAMETERS 

col 

33 

FIRST 

PROBE  TO  BE  EDITED 

(ID 

34 

SECOND 

i PROBE  TO  BE  EDITED 

(ID 

35 

THIRD 

PROVE  TO  BE  EDITED 

(ID 

36-37 

FIRST 

CHANNEL  TO  BE  EDITED 

(12) 

38-39 

SECOND 

| CHANNEL  TO  BE  EDITED 

(12) 

40-41 

THIRD 

CHANNEL  TO  BE  EDITED 

(12) 

64-65  FIFTEENTH  CHANNEL  TO  BE  EDITED  (12) 

THE  NED  OPTION  PERMITS  SELECTIVE  ZEROING  OF  DATA  BY  PROBE 
AND  CHANNEL  NUMBER 
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2.10.4  Output  Details 

The  output  listing  from  HIAC1D  is  shown  in  figure  14. 
The  first  line  of  output  shows  the  flight  date  (22  MAR  77) 
and  the  input  tape  (LTA027) . The  second  line  contains  the 
start  and  stop  times  over  which  this  data  had  been  averaged 
(15:  38:30-15:  38:39)  . 


TOTAL  LWC 
TOTAL  Z 
TOTAL  DO 
TOTAL  UN 
TOTAL  mK 


(liquid  water  content 
(radar  reflectivity 
(median  volume  diameter 
(number  density 
(LWC/ /z* 1000. ) 


gm/M*  * 3 ) 
mm**6/m**3) 
microns) 
NO/m**  3 ) 


are  all  calculated  from  the  cloud  and  precip  probes  only,  the 
scatter  probe  is  not  used  in  these  totals.  SF,  RF , ST  and 
RT  are  the  TWCI  (total  water  content  indicator)  components; 


PRES 

(pressure  millibars) 

ALT 

(altitude  kilometers) 

TAS 

(true  airspeed  meters 

per  second) 

TEMP 

(true  temperature  centigrade) 

DEWP 

(dewpoint  centigrade) 

JW 

(liquid  water  content 

gm/m**3) 

are  all  calibrated  VCO  values.  The  remainder  of  the  average 
summary  contains  the  particle  type  and  indicator.  The  indi- 
cator states  whether  the  particle  type  is  "fixed"  or  "changes". 
If  the  particle  type  had  changed  within  this  averaging  interval 
the  values  given  should  be  suspect.  The  "PRECIP  CHANNEL  ONE:" 
phrase  was  explained  earlier  in  this  section.  Total  counts 
is  the  sum  of  all  the  particles  counted  within  the  forty  five 
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channels  all  three  probes) . 

The  next  part  of  the  output  is  a channel  by  channel 
listing  of  the  average  interval.  The  column  under  "CL"  is 
the  channel  number.  The  columns  under  "DIAMETER"  are  the 
center  diameters  (in  microns)  of  any  particular  channel  during 
this  interval.  Counts  are  the  sum  of  all  the  particles  regis- 
tered in  each  channel.  The  columns  under  "NM-DENSITY"  contain 
the  normalized  number  density  of  each  channel.  Number  density 
(number  per  cubic  meter)  is  normalized  by  dividing  the  density 
by  its  associated  bar  width  in  mm.  The  columns  beneath  “LWC" 
contain  the  liquid  water  content  of  each  channel.  "NM-DENSITY" 
and  "LWC"  are  averages  over  the  whole  interval.  There  are 
three  sets  of  these  columns,  the  one  on  the  left  is  for  the 
scatter  probe,  the  one  on  the  right  for  the  precip  probe  and 
the  one  in  the  middle  is  data  derived  from  the  cloud  probe. 

Beneath  the  "COUNTS",  "NM-DENSITY"  and  "LWC"  columns  are 
the  totals  for  each  probe.  Lastly  beneath  these  probe  totals 
are  the  average  values  for  each  probe  of  mK,  Dq,  UN  and  Z. 
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2.11  LOOP 

2.11.1  Program  Description 

Program  LOOP  uses  the  input  data  as  catalogued  by  BOBIN 
and  calculates  MASS  and  MEDIAN  VOLUME  DIAMETER  for  a measured 
radar  reflectivity.  It  uses  empirical  distributions  of  dimen- 
sionless number  density  ND  4/M.  The  mathematical  description 
of  this  method  may  be  found  in  DPSI's  Final  Report  #TR-76-0182. 

This  program  is  designed  to  be  run  on  an  INTERCOM  terminal 
with  its  output  file  printed  at  a high  speed  terminal. 
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2.11.2  Operating  Instructions 

LOGIN , ID , PASSWORD , TTYNUMBER, SUP1 
ATTACH, LGO,LOOPX4 34 0 , ID=BELSKY 1 
LGOi 


The  terminal  interaction  is  fairly  self  explanatory  see 
figure  15.  Note  that  to  end  execution  a response  of  0.0  must 
be  entered  for  ENTER  RADAR  Z? 


2.11.3  Output  Details 

Figure  15  contains  a complete  sampie  session  showing  the 
typical  usage. 


Figure  15:  LOOP  Sample  Session 


LOOP 

PASS*?  1 

FLIGHT  DATE  C 50  CHARS.)*?  TESTING  DATA  FROM  E77-51  29  OCT  77 

i 

I 

MAXIMIW  D/D0  IN  INTEGRATION-  6.132 
ENTER  INSTRUMENT  MAXIMUM?  6.0 

CALCULATE  <C>  OR  GIVEN  CG)  MK?  G j 

MK*  ? 1 2.  0 

ENTER  RADAR  Z?  100.0 

ALPHA*  . 988483 1 BETA*  .9845171  CMK>=  12.000 

M <<MK)*CZR)**» 51  - 120.  000  ORIGINAL  D0  * .565 

OUTPUT  D0*  . 553 

OUTPUT  M-  1 20. 09  5 

OUTPUT  Z*  100.  001 

<MK)  AT  INFINITY*  12.010  I 

MEDIAN-  .551 
OBJECT  M-  1 19.  999 

OBJECT  Z-  99.998 

MAX.  DI  AM.  — 3.316  I 


ENTER  RADAR  Z?  0.  0 

TYPE  FOR  COMPUTING  CENTERS  BATCH*  TAPE6*  PRINT#  OKAP 
TYPE  FOR  LY  CENTERCNEXT  DOOR)  BATCH#  TAPE6#  PRINT*  AC*  OKAP 


END  LOOP 

2.768  CP  SECONDS  EXECUTION  TIME 
COMMAND-  BATCH#  TMPE6*  PRINT#  AC#  OKAP 
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2.12  BOBIN 

2.12.1  Program  Description  ■ 

Program  BOBIN  is  a data  input  program  that  prepares  data 
for  program  LOOP.  It  accepts  or  modifies  the  21  data  points 
that  specify  number  density  versus  size  (as  output  from  pro- 
gram KNPLT1D) . 

There  are  three  basic  functions  of  BOBIN.  Function  1 is 
to  input  a new  set  of  21  points.  Function  2 converts  these 
21  points  to  their  respective  N^  values  (number  density  in 
this  size  class)  and  then  permits  reentry  of  the  N^  values. 
Function  3 is  for  data  smoothing.  The  input  data  (21  points) 
is  smoothed  or  recalculated.  It  is  designed  to  be  run  on  an 
INTERCOM  terminal. 


126 


2.12.2  Operating  Instructions 

LOGIN , ID , PASSWORD , TTYNUMBER , SUP! 
REQUEST, NFILE , *PF! 

ATTACH , DATA, BOB I NX 4 340 , ID=BELSKY! 
DATA! 

INPUT  RESPONSES 
UNLOAD, FILE , NFILE! 

LOGOUT! 


2.] 2. 3 Output  Details 


A complete  output  sample  is  shown  on  the  following 
page  in  figure  16. 


« 
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* COMMAND-  ATTACH#  DAT  A,  BO  BI  NX  43  4 0,  I D=  BEL  SKY#  MR*  I 
PF  CYCLE  NO.  - 002 

COMMAND-  ATTACH#  LOOP#  LOOPX4340#  I D=  BEL  SKY#  M R=  1 

PF  CYCLE  NO.  ° 001 
COMMAND-  RECUEST#  N FILE#  *PF 

COMMAND-  DATA 

ENTERING  NEW  DATA  FILE?  NO 

PRINT  DATA  FILE?  YES 


X Y 


1 

.075 

3.  56218  170 

2 

. 150 

2.  70169  370 

3 

.250 

1.  63594480 

4 

.400 

.9  12788  68 

5 

. 600 

. 64542632 

6 

.800 

. 40244435 

7 

1.  000 

. 10986567 

8 

1.  200 

-. 18266973 

9 

1.400 

-.47527677 

10 

1.  600 

-.76783842 

11 

1.800 

- 1.  069  58020 

12 

2.  000 

- 1.  38387440 

13 

2.  200 

- 1.  70053140 

14 

2.400 

-2.  017  16370 

15 

2.  600 

-2.  33383220 

16 

2.800 

-2. 65045310 

17 

3.  000 

-2.9  67  1329  0 

18 

3.  200 

-3.  28  374250 

19 

3.400 

-3.  60043360 

20 

3.  600 

-3.9  1703180 

21 

3.800 

-4. 23373430 

22 

4.  000 

-4. 55043470 

Figure  16:  BOBIN  Sample  Output 


A B 

0.  00000000  0.  00000000 


CH/WGING  FITTING  FUNCTION?  NO 

CHANGING  A#  B COEFFICIENTS?  NO 

DO  YOU  WISH  TO  ENTER  NEW  N VALUES?  NO 

DO  YOU  WANT  THE  FITTING  FUNCTION  RECALCULATED?  NO 

NEW  DATA  FILE  COMFLETE 
rr  ? tv-  on  ew  ncxi-v/nAn 
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2.13  FLTPMS 

2.13.1  Program  Description 

FLTPMS  was  written  to  produce  a replacement  for  the 
PMS-1D  tape  (should  one  be  needed) . It  does  so  by  reading 
an  RTX/8  TU-10  flight  tape,  reformatting  the  data  to  be 
compatible  with  KN0LL1D,  and  writing  the  reformatted  data 
onto  a new  tape.  It  is  commonly  run  in  batch  mode. 

FLTPMS  requires  an  RTX/8  TU-10  flight  tape  as  its  input 
tape.  The  format  for  this  is  shown  in  Appendix  26. 

The  output  tape  produced  has  the  same  format  as  a PMS-lD 
data  tape  (see  Appendix  2)  . 
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2.13.2  Control  Cards 

JOBNM , TP 2 , CM6 5000, TlOO.  PROB  NO.  NAME 

PAUSE.  PLS  MOUNT  DISK  LYCPFI 
MOUNT , SN=LYCPFI , VSN=LYCPF I . 

SETNAME (LYCPFI) 

REQUEST, TAPE 1,S, HI, MT, RING, VSN=TAPEN01.  (FOR  KNOLL1D) 

REQUEST, TAPE 3, MT,S,VSN=TAPENO 3.  (FROM  RTX/8) 

ATTACH , A1 , FLTPMS , ID=FRANCIS 
FTN , I=Al , PL=999999 , OPT=2 . 

FILE (TAPEl , RT=U , BT=K , MRL=1024 , MBL=1024 , RB=1 , BFS=105 ) 

FILE (TAPE3 , RT=U , BT=K , MRL= 113  5 , MBL=1135 , RB=1 , BFS=116 ) 

LDSET , FILES=TAPEl/TAPE3 , PRESET=ZERO . 

LGO . 

7/8/9 

6/7/8/9 


2.13.3  Input  Cards 

none  required 
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2.14  EWERCOMP 

2.14.1  Program  Description 

Program  EWERCOMP  takes  the  EWER  VCO  values  from  the 
PMS-1D  buffer  and  calculates  a least  square  regression  line 
against  the  KNOLL ID  derived  total  liquid  water  content.  The 
output  consists  of  a scatter  plot  with  the  regression  line 
plotted  through  the  data,  and  an  output  listing  of  the  least 
square  coefficients,  RMS,  and  calculated  values  with  their 
variance  from  the  input  data. 
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2.14.2  Operating  Instructions 

LOGIN , ID , PASSWORD , TTYNUMB , SUP 1 
ATTACH , CRT , OFFLINECRT . 4- 
LIBRARY , CRT . i 
REQUEST, TAPE 39, *Q. 1 
DISPOSE , TAPE39 , *FL . + 

ATTACH , TAPE1 , PLOTTAPE , ID=GLASS . 1 (KN0LL1D  TAPE 2 ) 

ATTACH , LGO , EWE  RCOMP B I NX  2 876, ID=GLASS . + 

LGO.  4 
LOGOUT 
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2.15  EWERTIME 
.15.1  Program  Description 

Program  EWERTIME  was  written  to  plot  PMS-lD  total  liquid 
water  content  for  the  same  time  period  and  on  the  same  scale 
as  that  used  on  the  Aerospace  derived  plots.  This  permits  a 
visual  correlation  of  the  data.  Input  to  EWERTIME  is  the 
KNOLLlD  output  TAPE2 . Output  is  a CALCOMP  pen  plot.  EWERTIME 
can  be  executed  from  either  batch  or  INTERCOM. 
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2.15.2  Operating  Instructions 

JOBNM, CM60000 , T20 . ID 

ATTACH , PEN , ONLINEPEN . 

LIBRARY (PEN) . 

REQUEST, PLOT , *Q. 

DISPOSE, PLOT, *PL. 

PAUSE.  PLS  MOUNT  DISKLYCPFI 
MOUNT  SN=LYCPFI , VSN=LYCPFI . 

SETNAME (LYCPFI ) 

ATTACH , TAPEl , PLOTTAPE , ID=GLASS . (FROM 

ATTACH , LGO , EWERTIMEBIN , ID=GLASS . 

LGO . 

6/7/8/9 


NAME 


KNOLL1D ) 
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2 . 16  HIGHDEN 

2.16.1  Program  Description 

HIGHDEN  was  written  to  produce  scatter  plots  of  height 
(in  feet)  vs.  the  log  of  the  number  density  (N/M3)  measured 
during  a user  specified  interval  of  time.  It  can  prepare 
plots  from  either  a standard  INTERCOM  terminal  (by  disposing 
the  output  to  a plot  queue)  or  in  real  time  from  the  Tektronix 
graphic  terminal. 

HIGHDEN  uses  the  standard  KNOLLlD  plot  tape  which  has 
previously  been  cataloged  onto  the  disk  as  its  input  tape. 
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2.16.2  Operating  Instructions 


To  use  HIGHDEN  from  a standard  INTERCOM  terminal,  the 
user  must  execute  the  following  commands: 


LOGIN , NAME , PASSWORD , TTYNOI 
M,  PLS  MOUNT  DISK  LYCPFI4- 
MOUNT , SN=LYCPFI , VSN=LYCPFI I 
SETNAME (LYCPFI ) 1 

ATTACH , LIB , SMALLONLINE , SN=SHARED+ 
LIBRARY (LIB) V 
REQUEST, SMPLOT, *Q1 

ATTACH , TAPEl , KNOLL1DPLOTTAPE , ID=NAME I 
ATTACH , LGO , HIGHDENBIN , ID=LEACH1 
LGOl 

[INPUT  REQUESTS] 

DISPOSE , SMPLOT , SP 1 


To  use  HIGHDEN  from  the  Tektronix  plotting  terminal,  the 
user  must  execute  the  following  commands: 

LOGIN , NAME , PASSWORD , TTYNO 
M,  PLS  MOUNT  DISK  LYCPFI 
MOUNT, SN=LYCPFI , VSN=LYCPFI 
SETNAME (LYCPFI) 

ATTACH, LIB, TEKL IB, SN= SHARED 
LIBRARY (LIB) 

ATTACH , TAPEl , KNOLL 1DPLOTTAPE , ID=NAME 

ATTACH , LGO , HIGHDENBIN , ID=LEACH 

LGO 

[INPUT  REQUESTS] 
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2.16.3  Plotting  v , tions 

Regardless  of  which  type  of  terminal  is  being  used,  the 
program  will  begin  asking  questions  about  the  grapn  to  be 
plotted  after  the  LGO  command  is  executed.  The  user  must 
input  the  following  parameters: 

(1)  Whether  a standard  INTERCOM  terminal  or 
the  Tektronix  plotting  terminal  is  being 
used 

(2)  Which  probe's  number  density  (or  their 
total)  is  to  be  graphed 

(3)  Whether  normalized  or  unnormalized  density 
is  to  be  used 

(4)  The  date  of  the  flight 

(5)  Whether  the  time  is  to  be  taken  from  the 
aircraft  clock  or  the  PMS  clock 

(6)  The  starting  time  of  the  interval  of 
interest 

(7)  The  ending  time  of  the  interval  of  interest 

For  example,  the  terminal  session  in  Figure  17  will  pro- 
duce a scatter  plot  of  the  height  vs.  the  log  of  the  total 
unnormalized  number  density  measured  from  18:18:10  until 
19:58:10.  The  time  will  be  taken  from  the  PMS  clock.  The 
date  of  the  flight  was  February  26,  1978  and  the  plot  was 
being  prepared  from  the  Tektronix  plotting  terminal.  Figure  18 
shows  the  graph  produced  as  a result  of  this  session. 


LGO 

ENTER  TVPE  OF  TERMINAL 
(STANDARD  INTERCOM- 1 , TEKTRONIX-2  ) 
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Fiyure  17:  JIIGIIDKN  Sample  Session 
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3.1  KN2UTIL 

3.1.1  Program  Description 

This  program  provides  a means  of  verifying  the  correct 
operation  of  the  PMS-2D  devices.  It  also  produces  a data 
summary  that  is  useful  in  the  manual  particle  typing  required 
for  other  analyses.  KN2UTIL  reads  the  standard  nine  track 
PMS-2D  particle  image  tape  (Appendix  13  & 14).  The  program 
must  be  run  through  the  6600  batch  processor.  Output  is  in 
the  form  of  a line  printer  listing  and  either  CRT  or  line 
printer  generated  images. 

To  summarize,  the  following  results  are  produced  by  this 
program: 

1.  A tape  summary  of  all  data  recorded  that  includes 

a.  record  type  (slow  or  fast) 

b.  record  length 

c.  record  number  (absolute  and  by  type) 

d.  record  time 

2.  A data  listing  by  record  type  as  specified 

a.  slow  data  - line  printer 

b.  fast  data  (on  selected  device) 

1.  line  printer 

2.  35  mm  film 
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.1.2  Control  Ca  Is 

JOBNM, CM75000 , T400 , NTl . PROB  NO.  NAME 

PAUSE.  PLS  MOUNT  DISK  LYCPFI 
MOUNT , SN=LYCPFI , VSN=LYCPFI . 

ATTACH, CRT ,CRTP LOTS. 

LIBRARY, CRT.  1 

REQUEST, TAPE 3 9, *Q. 

DISPOSE, TAPE39,*FM. 

SETNAME (LYCPFI) 

VSN , TAPEl=TAPENO/NT . (TAPENO  IS  PMS-2D  DATA  ACQUISITION  TAPE) 

REQUEST , TAPE1 , PE , L , NR , NT . 

FILE (TAPE1 , RT=U , BT=K , MRL=  5576 , MBL=55  7 6 , RB=1 , BFS=56  0 ) 

ATTACH , LGO , KN2UTILBIN , ID=GLASS , MR=1 . 

LDSET , FILES=TAPEl , PRESET= ZERO . 

LGO. 

REWIND , SUM, SOUT , LOUT . 

COPY, SUM. 2 
COPY, SOUT. 3 
COPY, LOUT. 4 
EXIT (S) 

COPY, SUM. 2 
COPY, SOUT. 3 
COPY, LOUT. 4 
7/8/9 
DATA 
6/7/8/9 

1 REMOVE  WHEN  CRT  NOT  DESIRED  (ID  CARD  COL  1-3  MUST  BE  BLANK) 

2 REMOVE  WHEN  SUMMARY  NOT  WANTED 

3 REMOVE  WHEN  SHORT  RECORD  LISTING  NOT  WANTED 

REMOVE  WHEN  LONG  RECORD  LISTING  NOT  WANTED 
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3.1.3  Data  Cards 

CARD  1 ID  CARD 

cc  1-3  PEN  FOR  CRT  OUTPUT 

BLANK  FOR  LONG  RECORD  LISTING  ON  PAPER 

CC  11-16  TAPE  NUMBER 

cc  21-26  FLIGHT  DATE  DDMONYR  (NO  SPACES) 

CARD  2 OPTION  CARD 

cc  2-3  NUMBER  OF  END  OF  FILES  TO  PROCESS  (DEFAULT=1) 
cc  6-10  NUMBER  OF  ABSOLUTE  RECORDS  TO  READ  (DEFAULT=999999) 

CARD  3 SS  NAMELIST  CARD* 

cc  2-4  $ SS 

VALID  VARIABLES 

ARECB  ABSOLUTE  BEGINNING  RECORD 

ARECE  ABSOLUTE  ENDING  RECORD 

SRECB  BEGINNING  SLOW  RECORD  TO  BE  LISTED 
SRECE  ENDING  SLOW  RECORD  TO  BE  LISTED 
LRECB  BEGINNING  LONG  RECORD  TO  BE  LISTED 
LRECE  ENDING  LONG  RECORD  TO  BE  LISTED 

CARD  4 

CC  2-5  $END 

CARD  5 TIMEFLAG  NAMELIST  CARD** 

cc  2-9  $ TIME LAG 

VALID  VARIABLES 

TF  CHANGE  DEFAULT  TIME  FLAG  CODE 

(DEFAULT  IS  TF  = 0,1, 0,1, 0,1, 0,1) 


CARD  6 


cc 


2-5  $END 
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$ SS  cards  are  in  standard  NAMELIST  format.  When  the 
ARECB  and  ARECE  variables  are  specified  all  records 
whether  they  be  long  or  slow  that  are  between  records 
ARECB  and  ARECE  inclusively  will  be  listed.  (The  ARECB 
or  ARECE  variable  specifies  the  actual  tape  record  num- 
ber) . 

Specific  short  and/or  long  records  may  be  dumped  selec- 
tively by  using  the  SRECB,  SRECE,  LRECB  and  LRECE  variables. 
The  example  in  figure  19  depicts  the  correct  usage. 

If  the  CRT  option  was  selected  then  the  long  records 
will  be  put  onto  35  mm  film. 


$ TIMFLAG  cards  are  in  standard  namelist  format.  TF 
is  dimensioned  as  length  = 8.  Timing  word  is  indicated 
by  a bit  pattern  of  01010101  in  bits  1-8  of  any  scan, 
however  due  to  hardware  problems  it  may  be  something 
different.  The  TIMFLAG  cards  are  a means  of  changing 
the  flag  key  AFTER  it  is  known.  If  the  incorrect  pat- 
tern is  00101010  then  the  namelist  cards  should  read: 

@ $TIMEFLAG 

(§  TF=  0,0, 1,0, 1,0, 1,0, 

@ $END 


where  @ = blank  column 
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Record  type 

Record  number 
Slow  sequence 
Fast  sequence 


CASE  A:  To  dump  all  of  the  above  records  either  of  the 

following  methods  should  be  used. 

1.  records  1-20 

or 

2.  Slow  records  1-9 

and 

Fast  records  1-11 

USE:  @$SS  ARECB=1, ARECE=20  $END 

or 

@$SS  SRECB=1 , SRECE=9 , FRECB=1 , FRECE=9  $END 

CASE  B:  To  dump  the  records  indicated  by  the  aboe  rectangel 

use  either: 

1.  records  8-15 

or 

2.  Slow  records  5-7 

and 

Fast  records  4-8 
@$SS  ARECB=  8 , ARECE= 1 5 $END 
or 

@$SS  SRECB=5,SRECE=7,FRECB=4,FRECE=6  $END 
where  @ = blank  column 


Figure  19:  KN2UTIL  Sample  Input 
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3.1.4  Output  Derails 

Figure  20A  shows  the  tape  summary  that  is  produced  when- 
ever KN2UTIL  is  run.  The  listing  shows  the  count  of  each 
record  in  terms  of  both  total  records  and  record  type.  Each 
long  record  (indicated  by  an  L in  column  1)  has  a probe  iden- 
tifier (CL  for  cloud,  PR  for  precip)  shown  in  column  3.  All 
long  records  (particle  data)  must  be  547  words  in  length. 

The  short  records  (an  S in  column  2)  must  have  a length  of 
86  words.  A sync  byte  (should  be  1111),  a clock  time,  and 
an  elapsed  second  counter  are  all  shown  for  each  slow  record. 
Any  record  that  is  not  either  547  or  86  words  is  indicated 
by  an  X in  the  left  hand  column.  These  records  are  ignored 
by  all  other  PMS-2D  processing  programs. 

A short  record  summary  is  depicted  in  figure  20B.  These 
records  contain  ten  1-second  VCO  samples,  status  words  for 
each  probe,  and  five  percent  clock  rate  samples.  Also  there 
is  space  allocated  to  multiplex  15  PMS-lD  channel  counts  each 
second,  and  16  analog  data  channels  each  record.  At  the  pre- 
sent time  this  space  is  not  being  utilized.  The  eight  VCO 1 s 
that  appear  at  the  bottom,  are  calibrated,  and  averaged  (over 
10  seconds)  by  KN2UTIL;  they  appear  .for  convenience  only  and 
are  a means  of  visually  verifying  the  VCO  hardware. 

A long  record  line  printer  output  is  shown  in  figures 
20  C & D.  This  output  is  quite  lengthy,  at  one  line  per 
scan,  a three  record  set  uses  approximately  16  pages.  At 
the  beginning  of  each  record,  the  times  (clock  & elapsed) , 
record  count,  and  overload  status  are  shown.  The  elapsed 
time  (based  on  the  current  clock  sample  rate)  between  par- 
ticles is  also  shown.  Note  if  a record  appears  without  an 
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elapsed  second  count  between  each  particle,  there  is  a good 
chance  that  the  time  flag  may  be  incorrect  (see  Section  3.1.3 
for  a method  to  correct  this) . At  the  end  of  each  record 
(figure  20D)  the  elapsed  time  between  particles  is  summed 
and  printed. 

Figure  20E  also  depicts  the  long  records  using  a dif- 
ferent media.  Each  record  is  shown  on  one  35  mm  frame.  The 
information  is  identical  to  that  on  the  line  printer.  The 
big  advantage  to  this  media  is  the  amount  of  paper  used. 

Whereas  the  line  printer  output  has  a faster  (same  day,  usually) 
throughput  time . 
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3 . 2 TWODEE 

3.2.1  Program  Description 

This  program  is  concerned  with  transforming  the  32  diode 
scans  into  actual  particles  or  crystals.  The  technique  de- 
veloped uses  the  real  time  principle  of  a one-pass  calculation 
that  is,  while  the  particle  is  being  determined  certain  key 
parameters  are  also  being  calculated.  This  means  that  after 
the  initial  pattern  recognition  pass,  the  definition  of  the 
particle  is  fully  specified. 

Program  TWODEE  uses  the  standard  PMS-2D  particle  image 
tape  (Appendix  14  & 15 ) . A particle  tape(s)  (Appendix  13) 
is  produced  as  input  to  program  KN0LL2D.  An  output  listing 
that  contains  a summary  of  the  particle  tape  is  also  included. 
This  summary  can  optionally  be  selected  as  either  a short  or 
long  format.  The  short  format  includes  particle  header  infor- 
mation and  the  fundamental  parameters  of  each  particle. 
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3.2.2  Control  Cards 

JOBNM, CM65000 , T2000 , NTl . PROB  NO.  NAME 

ATTACH , LL , TIMELFT . 

FILE (TAPEl , RT=U , BT=K , MRL=5  576 , MBL=5  576 , RB=1 , BFS=560 ) 
PAUSE.  PLS  MOUNT  DISK  LYCPFI 
MOUNT, SN=LYCPFI ,VSN=LYCPFI . 

SETNAME (LYCPFI) 

ATTACH , LGO, TWODEEBIN, I D= GLASS , MR=1 . 

REQUEST , TAPE2 , *PF , SN=LYCPFI . 

VSN, TAPEl=TAPENO/NT. 

REQUEST, TAPEl, PE, L, NR, NT. 

LDSET , PRESET=ZERO , FILES=TAPE1 . 

LOAD , LGO , LL . * 

EXECUTE. 

CATALOG , TAPE 2 , TWODDATA , ID=GLASS . 

REWIND, TAPE 3. 

COPY, TAPE 3. 

EXIT (S) 

CATALOG , TAPE2 , TWODDATA , ID=GLASS . 

REWIND, TAPE 3. 

COPY, TAPE 3. 

7/8/9 

DATA  CARDS 
6/7/8/9 
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Data  Cards 

CARD  1 IDENTIFICATION 

CARD 

cc  1-6 

XYR-NN 

FLIGHT 

ID  (COL  1 MUST  BE 

E OR  L) 

cc  11-19 

DDMONYR 

FLIGHT 

DATE 

CARD  2 OPTION 

CARD 

cc  5 

NLEN 

0 = FULL  CALCULATION  OF 

MAXIMUM 

LENGTH 

1 = APPROXIMATION  OF  MAXIMUM  LENGTH 


CC 

6-10 

NEOF 

NUMBER  OF  END  OF  FILES 

TO  PROCESS 

cc 

11-15 

NOUT 

0 = SUMMARY  OUTPUT 

1 = FULL  OUTPUT 

cc 

16-20 

NRST 

NUMBER  OF  FEET  USED  IN 

PREVIOUS 

EXECUTION 

(USED  FOR  RESTARTING  JOB) 


CC  21-25  NTSW*  0 = USE  REAL  TIME  CLOCK 

1 = INPUT  TIME  OF  FIRST  SLOW  RECORD 

CC  76-80  NFEET  NUMBER  OF  FEET  ON  CURRENT  REEL 

(DEFAULT  2400) 

CARD  3-  N PASS  CARDS 

CC  2-9  HH : MM: SS  SAMPLING  START  TIME 

CC  12-19  HH:MM: SS  SAMPLING  STOP  TIME 

cc  25  QTYPE  = I IF  SNOW  OR  ICE  IN  PASS 

= R IS  ONLY  RAIN  IN  PASS 

* IF  NTSW=1  NEXT  CARD  MUST  BE  A TIME  CARD  IN  FORMAT 
@HH : MM : SS 


where  @ = blank  column 
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3.2.4  Output  De.  ils 

Figure  21A  depicts  the  short  particle  tape  summary  list- 
ing. There  is  one  row  of  information  per  record.  This  in- 
cludes: record  number  and  length,  clock  time,  record  number 
from  the  nine  track  tape,  slow  record  number,  sample  time  and 
rate,  aircraft,  probe  and  overload  information,  and  also  six 
VCO  counts. 

The  full  particle  tape  summary  (figure  21B)  includes  all 
of  the  above  in  addition  to  the  fundamental  parameters  for 
each  particle  contained  within  the  record.  In  this  listing 
the  header  information  is  shown  in  the  two  half-lines  at  the 
top.  Each  full  line  contains  the  particle  parameters. 


TWODEE/TWODDMP  Sample  Output 
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3 . 3 TWODDMP 

.3.1  Program  Description 

In  the  event  that  TWODEE  was  run  without  the  summary 
option,  program  TWODDMP  will  produce  a similar  summary 
without  actually  recreating  the  particle  tape.  This  saves 
considerable  computer  time.  Note  that  the  input  tape  to 
this  program  is  a particle  tape  (TWODEE  output)  and  not  the 
PMS-2D  acquisition  tape. 
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3.3.2  Control  Cards 

JOBNM, CM65000 , T100 , TPl . PROB  NO. 

PAUSE.  PLS  MOUNT  DISK  LYCPFI 
MOUNT , SN=LYCPFI , VSN= LYCPFI . 

SETNAME, LYCPFI. 

VSN , TAPE2=TAPENO . 

REQUEST, TAPE2  . 

ATTACH  LOG , TWODDMPBIN, ID-GLASS ,MR=1 . 

LDSET , PRESET= ZERO . 

LGO . 

7/8/9 
DATA 
6/ 7/ 8/9 


NAME 


i 
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3.3.3  Data  Cards 

Only  one  card  is  required.  A blank  in  column  1 will 
produce  the  long  form  output;  any  other  character  in  column 
1 may  be  used  for  the  short  format. 


3.3.4  Output  Details 

The  summary  output  from  this  program  is  identical  to 
the  summary  produced  by  program  TWODEE.  Refer  to  section 

3.2.4  for  a description. 
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3.4  KNOLL2D 

3.4.1  Program  Description 

Although  KNOLL2D  is  in  a developmental  stage,  and  is 
still  not  quite  considered  a production  program,  documenta- 
tion is  included  here  for  purposes  of  report  completeness. 
Eventually  it  will  replace  KNOLL1D  as  the  primary  post  pro- 
cessing program. 

The  program  utilizes  many  of  the  routines  and  algorithms 
of  KNOLL1D,  however  the  big  advantage  is  the  greatly  enhanced 
data.  Input  is  from  the  particle  tape  produced  as  output  from 
program  TWODEE.  The  output  includes  all  of  that  produced  by 
the  ID  system  plus  additional  tables  and  summaries.  It  is 
run  through  the  6600  batch  processor. 
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3.4.2  Control  Cards 

JOBNM , CM65000 , T4 00 . PROB  NO. 

PAUSE.  PLS  MOUNT  DISK  LYCPFI 
MOUNT , SN=LYCPFI , VSN=LYCPFI . 

SETNAME (LYCPFI) 

ATTACH , TAPE2 , TWODDATA , ID=GLASS , MR= 1 . 

ATTACH  f LGO , KNOLL2DBIN , ID=GLASS , MR=1 . 

LGO . 

7/8/9 

DATA  CARDS 
6/7/8/9 


NAME 
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. 3 Data  Cards 

CARD  1 FLIGHT  CARD 

cc  1-6  XYR-NN 

cc  11-19  DDMONYR 

CARD  2 OPTION  CARD 


CARD  3-N 

PASS  CARDS 

cc  2-9 

HH:MM: SS 

cc  12-19 

HH : MM: SS 

cc  24-25 

cc  29-30 

FLIGHT  IDENTIFICATION 
FLIGHT  DATE 

(SECONDS) 

KISTLER  PRESSURE 
= 1 BACKUP  PRESSURE 

START  OF  SAMPLING  TIME 
STOP  OF  SAMPLING  TIME 
PARTICLE  TYPE  FOR  CLOUD  PROBE 
PARTICLE  TYPE  FOR  PRECIP  PROBE 


cc  AVERAGING  INTERVAL 

cc  10  IPRS  = o 
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3.4.4  Output  Details 

The  first  type  of  output  produced  by  KNOLL2D  is  similar 

to  that  produced  by  KNOLLlD . All  references  made  are  to 

figure  22A. 

A)  The  number  of  one  second  data  samples  that  were  averaged 
to  make  this  table 

B)  the  start  and  stop  time  of  this  interval 

C)  flight  identification 

D)  particle  typing  indicators 

E)  the  channel  number  for  reference 

F)  the  center  diameter  for  this  channel  of  this  probe.  It 
is  calculated  by  melting  the  maximum  length. 

G)  the  normalized  density.  It  is  calculated  by  computing 
the  number  of  particles  that  would  be  detected  by  this 
channel  size  in  a cubic  meter  of  sample  volume.  Then 
for  comparison  with  the  other  channels  it  is  normalized 
by  dividing  it  by  the  channel  barwidth. 

H)  The  liquid  water  content  for  this  channel  above  are  only 
repeated  once  for  each  channel  of  each  probe.  F,  G,  and 
H,  cloud  on  the  left,  and  the  precip  probe  on  the  right. 

I)  The  set  of  calibrated  VCO  and  VCO  derived  values.  The 
basic  VCO  values  are  PRESSURE,  DEWPOINT  and  TRUE  AIRSPEED. 


* 


HEIGHT  is  calculated  from  PRESSURE.  TEMPERATURE  and  JW-LWC 
are  both  VCO 1 s that  are  adjusted  by  airspeed.  The  C AIRSPEED 
is  a calculated  alue  given  by  PRESSURE,  PRESSURE  GRADIENT 
and  TEMPERATURE. 

J)  The  column  under  REJ  is  the  number  of  particle  used  for 
computation  in  each  channel  when  end  rejection  is  applied. 

TOT  is  the  total  number  counted  in  this  channel  Thus 
TOT-RE J=  number  of  particle  touching  an  end  diode.  There 
is  one  table  for  each  probe. 

K)  Under  each  probe  is  a probe  summary  of  various  meteorological 
parameters.  There  is  also  a set  for  the  TOTAL,  this  is 

the  total  of  cloud  and  precip  combined 


1) 

M 

liquid  water  - 

content 

2) 

Z 

derived  radar 

reflectivity 

3) 

DO 

median  volume 

diameter 

4) 

MK 

ratio  of  M to 

the  square  root  of 

5) 

SAMPLE 

(SEC) 

how  many  seconds  elapsed  collecting  data  for 
this  average.  This  is  the  sum  of  the  timing 
marks  between  all  the  particles  that  made  up 
this  average. 
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There  are  five  different  ways  of  categorizing  two  dimen- 
sional particles,  by  maximum  length,  horizontal  feret  projec- 
tion, area,  average  projection  ratio,  and  equivalent  circle 
ration.  The  relationships  between  these  categories  is  il- 
lustrated by  KN0LL2D.  There  are  five  pages  of  distribution 
matrices,  two  per  page,  one  for  the  cloud  probe  and  one  for 
precip  (see  figures  22B  - 22F) . They  demonstrate  the  rela- 
tions between  any  two  of  these  categories.  These  matrices 
are  given  once  per  pass  and  are  therefore  pass  totals.  Re- 
ferences are  given  to  figure  22B. 


A)  The  pass  that  generated  this  data 

B)  Sampling  start  and  stop  times 

C)  The  x axis  parameter 

D)  which  probe 

l 

E)  particle  type  used 

F)  Total  number  of  sampling  seconds  (see  K-5  in  figure  22A) 

G)  Totals  of  this  row  and  which  "channels"  were  used  in  these 
totals,  values  are  sums  of  number  density 

H)  the  y axis  paramter 

I)  The  number  density  of  all  particles  fitting  in  this 

x,y  intercept.  This  number  is  given  in  scientific  notation, 
x.e.  uv  + z = uv.x  10  The  + sign  can  take  on  the  following 
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values  and  meaning: 

+ positive  exponent 

* positive  exponent  and  only  one  particle  made  up  this  entry 
- negative  exponent 

(underline)  negative  exponent  and  only  one  particle 
made  up  this  entry 

*J)  The  size  limits  of  this  row.  i.e.  channel  01  has  those 
particles  whose  maximum  length  is  one  or  two  diodes  long 

*K)  same  as  J for  the  columns 

L)  same  as  G for  the  columns 


* for  maximum  length,  horizontal  feret  projection,  and  area 
these  limits  are  number  of  diodes, 
for  the  ratios  they  are  non-dimensional  numbers 


The  last  part  of  a KNCLL2D  printout  is  a comparison  of 
LWC,  Z and  other  parameters  for  a complete  pass,  calculated 
by  different  methods.  The  method  used  for  the  first  section 
was  by  maximum  length.  Thus  the  first  of  these  three  compa- 
risons are  also  calculated  by  maximum  length,  the  second  by 
horizontal  Feret  projection,  the  third  and  last  is  calculated 
by  area.  The  output  is  shown  in  figures  22G  - 221  and  is 
identical  in  format  with  the  first  section  figure  22A.  The 
difference  between  the  outputs  is  that  these  pass  averages 
are  centered  one  per  page  and  have  the  plot  type  descriptor 
above  it. 
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Figure  221:  KN0LL2D  Sample  Output 
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3.5  SDIODE 

3.5.1  Program  Description 

Program  SDIODE  examines  the  data  for  hardware  anomalies 
such  as:  a diode  in  the  on  or  off  state  at  all  times,  an  in- 
correct time  flag  or  time  word.  Large  streakers  that  appear 
in  the  data  record  are  also  indicated.  The  technique  utilized 
is  similar  to  a longitudinal  redundancy  check  used  by  most 
tape  drive  controllers. 

Program  SDIODE  uses  as  its  input  the  9 track  PMS-2D  data 
acquisition  tape  as  written  by  the  on-board  PERTEC  recorder. 

The  operating  environment  for  SDIODE  is  from  card  input  through 
the  CDC  6600  batch  processor. 
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3.5.2  Control  Cards 

JOB , CM60000 , T400 , NT1 . ID  NAME 

PAUSE,  PLS  MOUNT  DISK  LYCPFI 
MOUNT , SN=LYCPFI , VSN=LYCPFI . 

SETNAME (LYCPFI) 

VSN , TAPE 1= TAP ENO/NT . 

REQUEST , TAPE 1 , PE , L , NR , NT . 

ATTACH , LGO , SDIODEBINX24  9 7 , ID=GLASS . 

FILE (TAPE1 , RT=U , BT=K , MRL=54  7 0 , MBL=54  70 , RB=1 , BFS=54  9 ) 
LDSET , FILES=TAPE2 , PRESET= ZERO . 

LGO. 


3.5.3  Input  Cards 

NONE  REQUIRED 


* 
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3.5.4  Output  Description 

One  line  of  output  appears  for  each  PMS-2D  tape  record. 
This  line  will  contain  32  sums  indicating  the  number  of  times, 
during  that  record,  a particular  diode  was  in  the  on  state. 

The  sums  for  diode  numbers  9 through  32  should  be  fairly 
random  with  values  between  200  and  500.  Excessively  high  or 
low  values  for  a particular  diode  indicate  a probable  error. 
Diodes  1 through  8 should  have  values  in  the  same  range,  how- 
ever, there  should  be  alternating  high/low  pattern  to  the 
eight  sums.  This  pattern  indicates  the  presence  of  a correct 
time  flag. 
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4 . 1 SPANDAR 

4.1.1  Program  Description 

Program  SPANDAR  takes  the  IBM  360  generated  radar  cor- 
relation tape  and  reformats  it  into  a CDC  6600  SCOPE-NOS/BE 
compatible  tape.  SPANDAR  can  only  be  run  under  the  CDC  6600 
batch  processor  via  card  input. 

Input  is  the  radar  data  tape  (Appendix  19)  containing 
DBZ  values  generated  during  an  aircraft  sampling  pass  in 
Wallops  Island  vicinity.  It  is  processed  by  the  Applied 
Physics  Lab  at  Johns  Hopkins  University  before  being  sent 
to  AFGL/LYC . 

There  are  five  different  output  options  available  from 
SPANDAR: 

1)  line  printer  plot  of  Z TAPE2 

2)  CDC  6600  data  tape  TAPE 3 

3)  data  summary  printout  TAPE7 

4)  tape  listing  (formatted)  TAPE8 

5)  punched  cards  TAPE 4 
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4.1.2  Control  Cards 

JOBNM, CM65000 , T50 , TP2 . * PROB  NO.  NAME 

REQUEST, TAPE 1 , L , MT , VSN=TAPENO.  (7  TRACK  - NO  RING) 

REQUEST, TAPE 3, MT, RING, VSN=TAPENO.*  (7  TRACK  - RING) 
ATTACH , LGO , SPANDARBIN , ID=GLASS , MR=1 . 

FILE (TAPEl , RT=U , BT=K , MRL=5  3 39 , MBL=53  39 , RB=1 , BFS=5  36 ) 
MAP, OFF 

LDSET , FILES=TAPE1 , P RES ET= ZERO . 

LGO. 


REWIND, TAPE8. 
COPY, TAPE 8. 
REWIND, TAPE2 . 
COPY , TAPE2 . 
REWIND, TAPE4 . 
COPY, TAPE 4. 
6/7/S/9 


FOR  TAPE  LISTING 


FOR  LINE  PRINTER  PLOT 


FOR  PUNCHED  CARDS 


* FOR  NO  OUTPUT  TAPE  REMOVE  REQUEST, TAPE 3, . . . CARD 
AND  CHANGE  TP 2 TO  TP1 


4.1.3  Data  Cards 

NONE  REQUIRED 


LI 
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4.1.4  Output  Description 

The  standard  SPANDAR  output  format  is  a summary  of  all 
the  data  on  the  converted  tape.  The  page  headings  specify 
which  aircraft  the  radar  was  tracking  and  each  block  under- 
neath specifies  one  correlated  radar-aircraft  track.  For 
instance,  in  figure  23A  the  C130  aircraft  was  tracked  on 
22  MAR  77  from  14:15:54  until  14:16:33  at  an  average  height 
of  1400  feet.  This  track  is  referenced  as  pass  3. 

Optional  TAPE 8 (figure  23B)  output  is  a complete  listing 
of  every  item  on  the  data  tape.  On  the  top  of  each  page  is 
an  expansion  of  the  pass  summary  block,  listing  aircraft  and 
times.  The  rest  of  the  page  is  devoted  to  the  tabular  values 
given  on  the  radar  tape: 

HH : MM: SS . F time  sample  was  collected  (once  a second) 

Z(DBZ)  radar  reflectivity 

EL (DEG)  radar  elevation  in  degrees 

AZ (DEG)  azimuth  in  degrees 

RSLRA(NM)  slant  range  to  the  aircraft  in  nautical  miles 

GRRA(NM)  ground  distance  to  the  aircraft  in  nautical  miles 

GR  RA (KM)  altitude  in  kilometer 

Optional  output  TAPE2  is  a line  printer  plot  figure  23C 
of  radar  reflectivity  (Z)  vs.  time.  The  range  of  Z is  -25.0 
DBZ  to  +25.0  DBZ  in  increments  of  .5  DBZ . The  program  divides 
the  data  into  one  minute  groups  with  a series  of  bars,  but 
does  not  cause  a break  in  the  graph.  The  time,  the  slant 
rante,  and  the  value  of  Z are  printed  to  the  side  of  each 
point  plotted. 
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Figure  23B:  SPANDAR  Sample  Output; 
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Figure  23C:  SPANDAR  Sample  Output 
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4 . 2 MORT 

4.2.1  Program  Description 

Program  MORT  reformats  the  radar  correlation  tape  pro- 
duced at  the  Kwajalein  Missile  Range.  This  tape  actually 
contains  data  from  two  different  radar  sources  (ALCOR  and 
TRADEX) . The  output  format  is  similar  to  that  produced  by 
program  SPANDAR.  This  permits  data  from  either  Wallops  or 
Kwajalein  to  be  used  interchangeably  by  program  RAPP.  Pro- 
gram MORT  must  be  run  through  the  CDC  6600  batch  processor. 
The  input  file  description  may  be  found  in  Appendix  22,  MORT 
input  format.  In  addition  to.  the  output  tape,  a complete 
tape  listing  is  also  produced.  The  output  file  description 
is  of  course,  identical  to  the  RAPP  input  tape  (Appendix  23) . 
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4.2.2  Control  Cards 

JOBNM,CM65000 , T100 , TP 2 . * ACT.  NAME 

REQUEST, TAPE 1 , MT , VSN=TAPENO . (7  TRACK  - NO  RING) 

REQUEST, TAPE3,MT, RING, VSN=TAPENO. * (7  TRACK  - RING) 

ATTACH  LGO , MORTBIN , ID=GLASS , MR=1 . 

MAP , OFF. 

LGO. 

7/8/9 

6/7/8/9 

* IF  AN  OUTPUT  TAPE  IS  NOT  NEEDED  CHANGE  TP 2 TO  TP1  AND 
REMOVE  REQUEST, TAPE 3, .. . CARD 


4.2.3  Data  Card 

NONE  REQUIRED 
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4 . 3 MORTDUMP/TTYMORT 
4.3.1  Program  Description 

As  a result  of  problems  trying  to  read  KREMS  tapes  from 
Kwajalein,  program  MORTDUMP  was  written  to  copy  these  radar 
correlation  tapes  to  disk.  Once  on  the  disk  the  INTERCOM 
text  editor  is  used  to  correct  the  errors  manually.  After 
the  errors  are  corrected  program  TTYMORT  is  used  to  format 
the  data  tape  as  if  it  had  been  processed  by  program  MORT. 
Only  the  operating  instructions  are  included  here,  for  a 
description  of  how  TTYMORT  operates  refer  to  program  MORT 
(Section  4.2).  TTYMORT  performs  the  same  function,  but  only 
from  the  terminal  under  INTERCOM. 
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4.3.2  Control  Cards 

JOBNM, CM65000 , T20 , TPl . PROB#  NAME 

REQUEST, TAPE 3, *PF. 

REQUEST, TAPEl,MT,VSN=TAPENO.  (7  TK  NORING) 

ATTACH  LGO , MORTDUMPBIN , ID=GLASS . 

LGO . 

CATALOG, TAPE 3 , KREMS , ID=FRANCIS 
6/7/S/9 


At  this  point  manually  edit  the  data  on  INTERCOM  by: 


ATTACH , TAPE , KREMS , ID=FRANCIS1 

REQUEST,  TAPE2,  *PF4- 

EDITOR+ 

E TAPE  S ! 

[EDIT  OUT  THE  ERRORS] 

S TAPEl  0 Nl 

ATTACH , LGO , TTYMORTBIN , ID=FRANCIS4- 
LGOi 

CATALOG , TAPE2 , RAPPTAPE , ID=FRANCISi 

AftSr  TAPE 2 has  been  catalogued,  it  is  now  ready  for 
standard  processing  with  program  RAPP. 
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4 . 4 RAPP 

4.4.1  Program  Description 

Program  RAPP  correlates  data  collected  from  the  PMS  ID 
and  2D  devices  aboard  the  C130-E  or  the  LEARJET  with  data 
acquired  from  a ground  based  radar.  The  correlated  data  is 
then  plotted  in  a variety  of  ways: 


1) 

aircraft 

"MK 

" vs.  radar  "Z" 

2) 

aircraft 

"Z" 

vs.  aircraft  "M" 

3) 

aircraft 

"Z" 

and  aircraft  "M" 

vs . time 

4) 

aircraft 

"Z" 

and  radar  "Z"  vs. 

time 

5) 

aircraft 

"MK 

" and  radar  "Z"  vs 

. time 

6) 

aircraft 

"Z" 

vs.  radar  "Z" 

7) 

aircraft 

"M" 

vs.  radar  "Z" 

M is  liquid  water  content  GM/M**3 
Z is  radar  reflectivity 
MK  is  defined  as  1000*M//Z 

There  are  two  tape  inputs  to  RAPP,  TAPE7  produced  by 
KN0LL1D  and  TAPE3  produced  by  either  SPANDAR,  MORT  or  TTYMORT. 
Since  RAPP  requires  tape  input  it  must  be  run  from  cards 
in  batch  mode. 

The  output  listing  consists  of  tabulated  data  points, 
correlation  percentages  for  shifts  of  aircraft-radar  match- 
ings of  +3  seconds,  and  least  square  regression  coefficients 
for  plots  1,2, 6, 7 listed  above.  The  graphical  output  uses 
35mm  film  generated  by  a CALCOMP  CRT  plotter. 
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4.4.2  Control  Cards 

JOBNM, CM67000 , T200 , TP2 . ACT.  NAME 

REQUEST, TAPEl,MT,VSN=TAPENO.  (7  TK-NORING) (FROM  KNOLL ID) 

REQUEST, TAPE3,MT,VSN=TAPENO.  (7  TK-NORING) (FROM  SPANDAR) 

ATTACH , LGO , RAPPBIN , ID=GLASS ,MR=1 . 

ATTACH , CRT , OFFLINECRT . 

LIBARARY (CRT) 

REQUEST, TAPE39, *Q. 

DISPOSE , TAPE39 , *FL . 

LDSET , PRESET= ZERO . 

LGO. 

7/8/9 

-DATA  CARDS- 
6/7/S/9 


AD-A064 


UNCLASSIFIED 


DIGITAL  PROGRAMMING  SERVICES  INC  WALTHAM  MASS  F/G  4/2 

DEVELOPMENT  and  APPLICATION  of  MATHEMATICAL  PROCEDURES  to  a VAR— ETC (U) 
71  LE  6ELSKY'  M V FRANCIS*  F B KAPLAN  F19628-76-C-0241 
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4.4.3  Data  Cards 


card  1 

CC  1-6  MISSION  ID  (XYR-NN) 

cc  8-13  A/C  TYPE  E.G.  C130-X 

CC  15-24  FLIGHT  DATE  DD  MON  YR  (LJ) 

cc  25  PROBE 

1 SCATTER 

2 CLOUD 

3 PRECIP 

4 TOTAL  (CLOUD+PRECIP ) 

CC  31-40  RADAR  OFFSET  (METERS)  F10 . 0 

cc  41-50  RADAR  CORRECTION  FI 0.1 

CC  55  KWAJ  ( 1-IF  KWAJALEIN  FLIGHT 

O-OTHERWISE)  - SEE  BELOW 

cc  56-65  DB1  (MINIMUM  RADAR  DETECTABLE  CONSTANT)  F10.0 

DBZ  MINIMUM  AT  1 KILOMETER 

CC  70  IQUIP  = 1 FOR  ID  DATA ( DEFAULT ) 

2 FOR  2D  DATA 

cards  2-N+l  (FOR  N PASSES) 

PASS  BEGINNING,  ENDING  CARDS 
FOR  EACH  DESIRED  PASS,  ENTER  CARD: 

cols  1-2  PASS  NUMBER 

cols  4-9  BEGIN  PASS  TIME,  AS  HHMMSS 

cols  11-16  END  PASS  TIME,  AS  HHMMSS 

NOTE:  TIMES  ARE  AIRCRAFT  TIMES  IF  IRADAR=0 , 

THEY  ARE  RADAR  TIMES  IF  IRADAR=1 

cols  17  I IF  ICE  IN  PASS;  W IF  WATER  IN  PASS 

cols  20-29  MINIMUM  M VALUE  (LWC)  (F10.0) 

cols  30-39  MAXIMUM  M VALUE  (LWC  IN  GRAMS/M**3)  (F10.0) 

col  45  RADAR  TYPE 

0 SPANDAR( DEFAULT) 

1 ALCOR  (KWAJ  MUST=1  ON  HEADER  CARD) 

2 TRADEX  (KWAJ  MUST=1  ON  HEADER  CARD) 


I 
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.4.4  RAPP  Output  Description 

The  first  page  of  a RAPP  output  (figure  24A)  gives  the 
aircraft  and  flight  identification  along  with  the  time  and 
date  of  processing. 

The  second  page  (figure  24B)  is  a pass  summary.  It 
specifies  the  pass  number,  actual  start  and  stop  time,  crystal 
type  for  the  pass  and  the  average  height  as  determined  by  the 
radar. 

The  first  analysis  section  (figure  24C)  is  the  aircraft 
derived  data.  The  probe  from  which  the  data  is  generated  is 
listed  on  the  first  line.  The  first  column  is  the  time  of 
the  sample.  The  columns  labelled  "AVE  Z"  and  "AVE  M"  are  the 
data  points  as  calculated  and  averaged  by  KN0LL1D.  The  columns 
labelled  "M" , "Z"  and  "MK"  are  the  values  interpolated  by 
Newton's  Fundamental  Formula.  They  represent  LWC,  radar  re- 
flectivity and  K factor. 

The  next  section  (figure  24D)  is  simply  a listing  of  the 
radar  tape  (Appendix  23) . The  additional  value  " DBZMIN"  is 
the  minimum  detectable  radar  signal  at  1 nautical  mile.  Any 
DBZ  value  below  this  value  will  be  rejected. 

The  last  two  output  pages  (figure  24E  & 24F)  are  the 
correlation  summary.  The  first  section  lists  what  data  cor- 
relation result  from  shifting  the  data  +3  seconds.  In  this 
example  the  best  correlation  is  no  shift.  As  each  plot  is 
made  the  following  parameters  are  printed. 
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1)  the  plot  being  made 

2)  the  number  of  prints  curve  was  fit  to 

3)  the  x and  y log  mean  and  standard  deviation 

4)  the  regression  on  x and  the  on  y,  with  coefficients 

5)  the  regression  with  the  lowest  RMS  is  printed  as  the 
one  used 

6)  the  correlation  of  the  data 
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~ figure  25B:  RAPP  Sample  Plots 
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5.1  RTX/8 

5.1.1  Program  Description 

The  RTX/8  program,  which  runs  on  the  DEC  PDP-8/E,  is 
actually  an  operating  system,  rather  than  an  executing 
program.  It  has  built  into  it:  the  capability  to  take 
inputs  from  external  devices,  process  the  data,  record  it 
on  magnetic  tape,  print  data  on  the  printer,  and  plot  on 
the  Tektronix  screen,  all  within  the  constraints  of  real 
time . 


There  are  two  versions  of  RTX/8;  one  for  the  lab 
(Airpl  = 0)  and  one  for  the  C-130  (Airpl  = 1) . The  dif- 
ference is:  on  the  plane,  the  printing  device  is  the  GE 
terminet,  in  the  lab,  it  is  the  ODEC  line  printer. 

RTX/8  is  also  capable  of  time-sharing  user  programs. 

The  following  is  a list  of  RTX/8  user  programs. 

MNSI*  System  interpreter.  This  activates  the  keyboard, 

and  starts  and  stops  other  user  programs 
TAPE  Writes  the  contents  of  the  RTX/8  buffer  onto  magtape 

PRINT  Prints  selected  parameters  on  the  line  printer 

UNIT*  A sub-user  initiated  by  PRINT,  which  calibrates  the 

raw  VCO  data  and  sends  it  to  the  print  buffer 
LWCD  Computes  liquid  water  content  from  raw  probe  counts 

PLOT  Plots  selected  parameters  on  the  Tektronix  screen 


* SPECIAL  USER  PROGRAMS  (NOT  EXECUTABLE  VIA  KEYBOARD) 
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.1.2  Operating  Instructions 

1)  mount  system  tape  on  dec tape  0. 

2)  turn  Tektronix  screen  and  hard  copy  unit  on 

3)  key  7470  into  switch  register 

4)  Press: 

a.  ADDR  LOAD 

b.  EXT  ADDR  LOAD 
C • CLEAR 

d.  CONTINUE 

Computer  will  respond  with  a dot  on  the  screen 

5)  type:  R RTX8  (must  include  a space  between  R and  RTX8) 

6)  RTX8  will  clear  the  screen,  type  current  version  number 
and  ask  for  today's  date  in  the  form  DDMMYY 

7)  put  console  knob  on  AC.  This  will  display  a rotation 
of  lights  on  the  computer. 


Commands  recognized  by  RTX8 


CTRL/C:  Causes  to  be  printed  on  the  screen  and  get 

the  attention  of  RTX8  so  programs  may  be  turned 
ON  and  OFF 

CTRL/E:  Clear  and  reset  the  screen  (If  PLOT  is  running, 

a new  graph  will  be  started.  NO  HARD  COPY) 
CTRL/P:  Causes  hard  copy,  clear  and  reset  the  screen. 

(If  PLOT  is  running,  a new  graph  will  be  started) 
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5.2  User  Program  LWCD 
5.2.1  Program  Description 

LWCD  computes  the  total  liquid  water  content  from  the 
raw  probe  counts.  Additionally,  it  performs  several  auxil- 
iary calculations.  These  include:  radar  reflectivity  (Z)  , 
for  factor  (F) , MK  and  stability  factor  (S) . 'LWC' , 'Z', 

'F',  'MK'  and  'S'  are  sent,  in  ASCII,  to  the  output  buffer 
used  by  program  PRINT.  They  are  also  sent  fixed  point,  single 
precision  (with  a scaling  factor)  to  the  PLOT  program. 


I 
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5.2.2  Operating  Instructions 

TYPE:  ON , LWCD1 

VELOCITY: 

HARDWARE  INPUT  OKAY? 


If  the  true  airspeed  is  to  be  used,  type  "Y". 
board  is  desired,  type  "N".  If  "N" , the  system  re 


If  key- 
sponds  with 


VELOCITY: 

Type  the  velocity  in  knots 

PARTICLE : 

HARDWARE  INPUT  OKAY? 

If  "N",  RTX-8  types: 
PARTICLE  TYPE  (1-5) ? 


Select  the  proper  particle  type: 


1 = RAIN 

2 = WET  SNOW 

3 = LARGE  SNOW 

4 = SMALL  SNOW 

5 = BULLET-ROSETTES 


PROBE :C,P,  or  B? 


Select  cloud  probe  (C) , precipitation  probe  (P)  or 


both  (B) . 


; . - - - - 
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5 . 3 User  program  TAPE 
5.3.1  Program  Description 

Program  TAPE  assumes  the  VCO  input  in  the  form  shown  in 
Appendix  24  and  it  assumes  the  PMS  data  is  coming  in  at  a 
rate  of  256  characters  (4  bits  each)  per  second.  The  output 
on  magnetic  tape  is  in  the  form  shown  in  Appendix  26. 

TAPE  can  also  read  PDP-8  flight  tapes  in  the  Lab.  To 
do  this,  4 locations  must  be  manually  changed  using  the 
switches  on  the  front  panel  of  the  computer. 


LOCATION 

FROM 

TO 

10644 

7650 

7200 

10711 

0001 

0002 

05032 

1725 

7200 

05033 

3725 

7200 
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5.3.2  Operating  Instructions 

1)  Put  a tape  on  the  TU-10 

(Make  sure  tape  drive  is  on-line  and  tape  has 
write  ring) 

2)  TYPE: 

ON, TAPE i 

To  terminate  this  program 

TYPE: 


OFF, TAPE I 
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5.4  User  program  PLOT 
5.4.1  Program  Description 

This  user  program  plots  the  data  (see  figure  26)  output 
by  the  LWCD  program  on  the  Tektronix  CRT.  This  includes  LWC, 
F,  S,  MK  and  Z.  The  operator  can  choose  any  one  of  these 
parameters  to  be  plotted;  if  M is  chosen,  a choice  of  three 
ranges  is  available.  The  axes  are  drawn  on  the  right  half 
of  the  CRT,  leaving  the  left  side  for  operator  dialog.  Each 
plot  is  labelled  with  function,  time,  maximum  value  and  date 
so  it  can  easily  be  identified.  Each  tick  mark  on  the  abcissa 
represents  one  minute.  On  the  ordinate  a tick  represents 
one-fifth  of  the  maximum.  After  five  minutes,  the  plot  fills 
up  and  automatically  a hard  copy  is  made,  the  screen  erased 
and  r.ew  axes  are  drawn  and  labelled.  A new  plot  and  hard 
copy  is  also  generated  whenever  the  screen  is  erased  by  op- 
erator dialog  filling  the  screen,  or  a control/P  function. 
Control/E  will  only  cuase  a new  plot  on  the  screen  without 
a copy. 

If  any  data  exceeds  the  plot  maximum,  an  interpolation 
is  done  and  a line  drawn  to  where  the  line  would  have  left 
the  plot  area.  This  also  occurs  when  the  data  returns  to 
the  range  of  the  plot. 


6JUN78 


LOG  <S> 

07-3O-21  TO  07-33-21 


Figure  26:  User  Program  PLOT  Sample  Plot 
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5.5  User  Program  PRINT 
5.5.1  Program  Description 

This  user  program  sends  one  line  of  information  to  the 
GE  Terminet  (or  ODEC  in  Lab)  printer  each  time  the  Knollenberg 
buffer  fills  (every  four  seconds) . This  contains  calibrated 
VCO  readings,  a time  code,  LWCD  results  and  a particle  dis- 
tribution. A brief  description  of  each  of  the  output  para- 
meters follows  (see  figure  27)  . 

PT  Particle  type  (R=rain,  W=wet  snow,  L=large  snow, 

S=small  snow,  B=bullet  rosettes 
ET  Elapsed  time  as  read  from  PMS-lD  buffer 

TIME  From  Stancil  Hoffman  time  code  generator.  Hours, 

minutes  and  seconds  (HH:MM:SS)  are  included  on 
line  1,  seconds  only  appear  on  lines  2-15 
ALT  Altitude  in  kilofeet:  calculated  as  a fifth  degree 

polynomial  from  the  Kistler  pressure  reading 
TEMP  Temperature  in  degrees  centigrade 

MAGH  Magnetic  heading  in  degrees  (0-136°N) 

DEWPOINT  in  degrees  Centigrade 

TAS  True  airspeed  as  read  from  the  TAS  computer 

JW-LWC  reading  from  the  Johnston-Williams  device 
EWER  LWC  reading  as  output  from  the  NOVA  computer 

LWC  liquid  water  content  re 

Z radar  reflectivity 

MK  M//Z  ratio 

F form  factor  computation 

S stability  factor 
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The  particle  distribution  is  displayed  as  one  minute  sums 
for  the  15  PMS  size  channels.  Note  that  the  45  sums  that  ap- 
pear are  sums  of  a one  minute  interval  ending  at  the  time 
associated  with  the  data.  That  is  in  the  output  format  shown 
in  figure  27,  the  first  set  of  sums  are  for  particle  during 
07:33:01  to  07: 34:00 
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5.5.2  Operating  Instructions 

1.  Turn  Terminet  printer  ON 

(The  switch  is  on  the  back  and  at  the  right  as  you  face 
it. ) 

2.  Press  the  LOCAL  button  hit  a carriage  return  and  make 
sure  there  is  a 1 in  the  print  column  light. 

3.  Position  paper  so  it  will  begin  printing  at  the  very  top 
of  the  page. 

4.  Press  the  ON  LINE  button,  the  READY  light  should  be  on 
if  everything  is  all  set. 

5.  To  activate  the  PRINT  program  type 

ON, PRINT 


NOTE:  The  PRINT  program  will  only  print  LWC,  Z,  MK,  F and  S 
if  the  LWCD  program  is  on.  Otherwise  only  the  VCO  and 
probe  counts  data  will  be  printed. 
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6.1  QWIK4 

6.1.1  Program  Description 

QWIK4  performs  a quick  look  type  dump  of  all  ID  PMS  and 
PDP8  data  tapes  from  the  C130E.  The  dump  printed  by  the  pro- 
gram includes  time  printouts,  raw  and  calibrated  VCO's,  total 
probe  counts,  and  certain  values  derived  from  this  data. 

The  program  may  be  run  in  two  different  modes.  The  first 
mode  dumps  the  entire  tape  printing  averages  over  a specified 
interval.  The  second  mode  allows  flight  time  or  elapsed  time 
to  be  used  to  locate  a record.  This  feature  will  be  useful 
in  locating  the  start  of  a sampling  run  or  any  time  interval 
of  interest. 


NOTE: 


There  are  two  versions  of  QWIK4 , one  for  the  plane 
and  one  for  the  Lab. 
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6.1.2  Operating  Instructions 

Mount  the  magtape  to  be  dumped  on  unit  0 with  the  write 
enable  ring  removed  to  prevent  any  possible  correuption  of 
data. 


Turn  on  the  GE  Terminet  printer  (or  ODEC  if  in  the  Lab) , 
press  'ON  LINE'  button  and  be  sure  that  the  paper  is  free  to 
feed. 


Respond  to  the  OS/8  dot  (.)  by  typing  "R  QWIK4 " (return). 

When  QWIK4  is  started  it  will  print 

0 = TTY,  1 = ODEC.  THEN  PRESS  CONTINUE. 

QUICK-LOOK  PROGRAM 

If  QWIK4  is  being  run  in  the  Lab,  it  will  type: 

"SET  SWITCH  #11. 

This  will  direct  output  to  DECwriter,  or  line  printer 

FLIGHT  TAPE  (F)  OR  KENNEDY  TAPE  (K) ? 

Respond  by  typing  F if  the  magtape  is  a PDP8  generated 
flight  tape  or  K if  the  tape  is  from  Kennedy  tape  recorder 
in  the  PMS  system.  Terminate  input  with  a return. 

NOTE:  all  inputs  to  QWIK4  are  terminated  with  a return. 


FLIGHT  NO. , DATE? 
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Enter  the  flight  number  and  date  or  a string  of  up  to 
30  characters  which  identifies  the  tape.  The  program  then 
asks 


DATE  PROCESSED , COMMENTS  ? 

Enter  date,  and  any  comments  (30  characters) 

ALL  (A)  OR  SPECIFIED  (S)  TIMES? 

Respond  with  A to  dump  the  whole  tape  starting  at  what- 
ever record  the  tape  head  is  currently  at.  The  program  will 
not  rewind  the  tape  to  the  load  point  before  dumping.  This 
is  useful  if  the  program  fails  to  locate  a specified  time 
interval.  The  operator  may  position  the  tape  to  approximately 
where  the  desired  time  would  be  using  the  off  line  controls 
on  the  TUlO . The  A option  may  then  be  used  to  determine  what 
was  recorded  on  the  tape.  There  is  one  small  problem  asso- 
ciated with  this  procedure.  When  the  off  line  controls  are 
used  to  forward  space  or  backspace  the  tape,  the  TUlO  does 
not  stop  on  an  even  record  boundary.  Therefore  the  first 
few  seconds  of  data  produced  by  the  program  should  be  ignored 
if  the  tape  has  been  moved  with  the  off-line  controls.  The 
program  will  detect  this  error  condition  and  print 

! ! ! INCORRECT  MAGTAPE  RECORD  LENGTH ! ! ! 

This  printout  is  normal  after  the  tape  has  been  moved  with 
the  off-line  controls,  however  it  should  only  print  once . 

The  first  read  done  by  the  TUlO  will  reposition  the  tape 
head  to  the  beginning  of  a record  and  subsequent  records 
should  be  readable  witr  generating  this  error.  If  this 
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printout  repeats  continuously  the  tape  is  not  in  the  correct 
format  and  any  results  printed  by  the  program  will  be  incorrect. 

The  S option  allows  the  operator  to  specify  a time  in- 
terval on  the  tape  to  be  dumped.  If  this  option  is  selected 
the  program  then  prints 

CLOCK  (C)  OR  ELAPSED  (E)  TIME? 

Type  C to  use  the  flight  time  clock  to  look  up  records  or  E 
to  use  the  PMS  buffer  elapsed  time.  The  program  will  attempt 
to  find  the  specified  time  interval  no  matter  where  it  is  on 
the  tape  and  will  minimize  tape  motion  as  much  as  possible. 

If  the  tape  has  gone  beyond  the  time  interval  desired  the 
program  will  backspace  as  necessary  to  rewind  the  tape  to 
the  load  point  before  each  run. 

After  the  operator  has  specified  which  clock  to  use  the 
program  will  print 

START  TIME 

Respond  by  typing  the  starting  time  for  the  dump  in  the 
form  SSSS  to  specify  a four  digit  elapsed  time  in  seconds  or 
HH : MM: SS  to  specify  a flight  time.  The  program  will  then 
print 

STOP  TIME 

Respond  with  the  time  to  stop  the  dump  at  in  the  same 
form  used  for  the  starting  time 


The  program  now  asks 


AVERAGING  INTERVAL  (SECS) 

This  question  is  printed  for  both  the  specified  time  interval 
and  dump  all  times  mode.  Respond  with  the  number  of  seconds 
to  be  represented  in  each  printout. 

The  program  then  begins  calculating  averages  over  the 
specified  time  interval  and  printing  a report  for  each  inter- 
val. Contained  in  each  printout  are  the  start  and  stop  flight 
and  elapsed  times  for  the  interval,  average  VCO  values,  average 
values  of  parameters  derived  from  the  VCO's,  and  total  probe 
counts . 

Some  sample  printout  is  contained  in  figure  28. 

The  program  may  be  restarted  from  the  beginning  by  set- 
ting the  switches  to  0200  and  pressing  HALT,  ADDR  LOAD,  EXT 
ADDR  LOAD,  CLEAR,  CONTINUE. 

Alternatively  the  program  can  be  restarted  from  the  "ALL 
(A)  OR  SPECIFIED  (S)  TIMES?"  question  by  setting  the  switches 
to  0000  and  pressing  HALT,  ADDR  LOAD,  EXT  ADDR  LOAD,  CLEAR, 
CONTINUE. 
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The  following  error  messages  may  print  out  during  execution 
of  QWIK4. 


printout 

possible  cause  (s) 

recommended  action (s) 

TIME  INTERVAL  OUT  OF 
RANGE 

time  specified  not 
on  tape 

use  "A"  option  to  see 
what  is  on  tape 

**END  OF  TAPE** 

time  specified  not 
on  tape 

rewind  tape  and  try 
another  time  interval 

! I ! MAGTAPE  PARITY 

ERROR! ! ! 

bad  tape,  bad  TU10 

try  another  tape; if 
the  printout  is  in- 
frequent it  may  be 
ignored; however  if 
data  looks  bad  this 
is  the  cause 

! ! ! INCORRECT  MAGTAPE 
RECORD  LENGTH! ! ! 

bad  tape  (unlikely) 
question  asking  if 
tape  was  Kennedy 
tape  (PMS ) or 

Flight  tape  (PDP8) 
was  answered  in- 
correctly tape  head 
was  not  on  even  re- 
cord boundary  due 
to  operator  moving 
it  using  off-line 
controls 

try  another  tape  re- 
start program  at  0200 
and  answer  question 
correctly 

no  action  necessary 
simply  ignore  results 
of  first  averaging 
interval 

Some  errors  will  not  be  detected  by  the  program.  Following 
is  a list  of  the  conditions  and  recommended  actions. 


disposition 

possible  causes (s) 

recommended  action (s) 

TU10  "rocks"  and  will 
not  locate  time  speci- 
fies 

bad  times  on  tape 

use  "A"  option  to 
dump  tape 

program  will  not  load 

bad  system  tape 
unknown 

try  backup  system  tape 
notify  DPSI 

program  "hangs  up" 

unknown 

record  AC,  PC,  MQ  and 
notify  DPSI 

program  halts 

unknown 

record  AC,  PC,  MQ  and 
notify  DPSI 

I'"  I M 


' I t-  M*,P  IlM 
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Figure  28:  QWIK4  Sample  Output 
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6 . 2 KNMON 

6.2.1  Program  Description 

Program  KNMON  was  written  for  testing  and  verification 
of  the  PMS-1D  interface.  It  also  allows  the  associated  M1703 
card  to  be  checked  out.  Any  desired  channel  selected  by  the 
operator  may  be  monitored.  Every  second  the  PMS-lD  system 
send  64  four  character  words  to  the  computer.  Actually  the 
PDP-8E  receives  its  data  one  character  at  a time  in  1/2S6 
second  intervals.  However  it  takes  one  second  for  the  64 
words  to  become  available.  Each  second  the  value  of  the  se- 
lected channel  is  printed.  Refer  to  Appendix  3 for  the 
table  of  PMS  channels. 


NOTE:  KNMON  can  only  function  correctly  when  run  on  the 
plane 
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6.2.2  Operating  Instructions 

1.  Respond  to  the  OS/8  dot  by  typing  "R  KNMON"  (return) 

2.  After  KNMON  is  started  it  will  print 

KN  CHANNEL 

Respond  by  typing  the  Knollenberg  channel  to  monitor 
(from  1 to  64)  followed  by  a return. 

3.  To  select  another  channel  set  SR  to  0200  and  press  'HALT', 
ADDR  LOAD,  ' EXTD  ADDR  LOAD',  'CLEAR',  'CONTINUE'  then 
continue  from  step  2. 
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6 . 3 PLOT 

6.3.1  Program  Description 

Program  PLOT  was  written  for  calibration  and  general 
plotting  at  LYC.  A modified  version  has  also  been  generated 
for  use  on  the  airplane  with  input  on  the  GE  keyboard  and 
output  on  the  CRT  and  GE  printer. 

Program  PLOT  will  ( 

(a)  plot  an  x-y  table  on  the  Tektronix  plotter;  the  table 

is  inputted  at  the  DECwriter  (or  GE  keyboard  if  on  plane) . • 

If  desired,  the  x and/or  y values  can  be  modified  by  log- 
ging them  in  order  to  produce  a linear-linear,  a linear- 
log,  a log-linear  or  a log-log  plot. 

(b)  generate  the  x values  automatically  after  a specification 
of  the  first  x and  x-step  is  given 

(c)  allow  the  user  to  select  the  low  and  high  values  of  x 
and  y to  be  used  on  the  plot 

(d)  plot  each  point  + ; the  user  can  choose  whether  or  not  to 
connect  the  plotted  points  with  a line 

(e)  label  the  plot  with  an  alphabetic  description 

(f)  plot  a least  square  best  fit  curve  to  the  data  (first  or 

second  degree) , and  print  out  the  fitting  function  J 


(g)  generate  a table  of  deviations  of  the  least  square  fit 
to  the  original  data,  and  calculate  the  RMS  error 
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(h)  allow  the  user  to  modify  the  data,  limits  and  descriptions 
and  replot,  with  a new  least  square  fit,  without  having 
to  retype  the  entire  x-y  table. 

An  example  of  the  plot  created  by  the  program  will  be 
found  in  figures  29  & 30. 
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6.3.2  Operating  Instructions 

In  order  to  execute  program  PLOT,  and  perform  the  many 
optional  capabilities,  the  following  step-by-step  procedure 
should  be  executed: 

1.  Mount  DECtape  containing  'PLOT'  on  unit  0,  write  enabled 
and  remote 

2.  Switches  set  to  7470 

3.  Press  ADDR-LOAD,  EXT-ADDR-LOAD,  CLEAR,  CONT 

4.  Computer  responds  with  a dot 

5.  Type  R PLOT 

6.  Within  22  seconds  the  computer  responds  with 

PLOTTING  PROGRAM 

7.  IS  X ON  LOG  SCALE? 

If  the  x-data  is  to  be  logged  base  10,  answer  Y otherwise  N 

8.  IS  Y ON  LOG  SCALE? 

If  the  y-data  is  to  be  logged  base  10,  answer  Y otherwise  N 

9.  IS  DELTA  X CONSTANT? 

If  the  x-data  (before  taking  logs)  is  equally  spaced, 
answer  Y,  otherwise  N 

10.  If  delta  x is  constant,  the  computer  will  request  X-START 
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and  DELTA-X,  as  X-START? 

Type  the  starting  value  and  return, 

DELTA-X? 

Type  the  difference  between  successive  x-values  and 
return 

11.  The  computer  responds  with 

ENTER  X , Y TABLE 

or  ENTER  Y TABLE  (for  constant  delta-x) 

12.  For  each  value  input  the  computer  will  first  print  the 
index  number  ( 1 , 2 , 3 , . . . etc . ) . If  only  y-values  are 
being  inputted,  the  user  will  type  the  y-value  followed 
by  a return;  if  both  x and  y values  are  being  inputted, 
the  user  will  type  the  x-value,  followed  by  a comma, 
followed  by  the  y-value  and  a return. 

13.  When  all  values  have  been  inputted,  type 

1E35  return 

14.  The  computer  will  respond  with 

LIMITS  OF  X ARE  xxxxx,  xxxxx 

which  indicates  the  low  and  high  values  of  the  x-table. 
Then  it  prints 

TYPE  LIMITS  TO  USE 

The  computer  expects  two  numbers,  separated  by  a comma, 
the  first  number  is  the  value  of  x at  the  left  end  of 
the  plot;  the  second  is  the  -value  of  x at  the  right  end 


m 
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of  the  plot.  There  will  automatically  be  10  divisions 
in  this  range.  If  the  first  limit  typed  is  unacceptable 
(larger  than  the  limit  found)  or  if  the  second  limit 
typed  (smaller  than  the  limit  found)  is  unacceptable, 
step  14  is  repeated. 

15.  Step  14  is  done  for  Y 

16.  The  computer  will  print 

CONNECTING  LINE? 

If  a line  is  desired  between  points  in  the  order  they 
were  typed  in,  answer  Y.  An  answer  of  N will  eliminate 
the  connecting  line. 

17.  The  computer  will  responde 

TURN  ON  PLOTTER,  TYPE  DESCRIPTION 

Type  a description  to  be  printed  on  the  bottom  of  the 
plot.  This  description  is  limited  to  67  characters. 

A runout  will  respond  with  a carriage  return,  and  the 
user  should  type  the  entire  description  over  again. 

Make  sure  plotter  is  on,  and  "on-line"  before  pressing 
return. 

18.  The  plot  will  be  made.  At  its  conclusion,  the  computer 
will  ask 

LEAST  SQUARE  FIT? 

If  a curve  fitting  the  data  is  desired  answer  Y;  other- 
wise answer  N and  proceed  to  step  22. 


233 


19.  The  computer  responds 

DEGREE  (1  OR  2)? 

Degree  1 is  a straight  line  fit;  degree  2 is  a quadratic 
fit.  Respond  either  1 or  2.  Any  other  response  will 
force  step  19  again. 

20.  The  fitting  curve  will  be  plotted.  The  least  square  poly- 
nomial will  be  printed 

(a)  at  the  top  of  the  plot,  and 

(b)  on  the  DECwriter  (or  GE  printer) 

21.  The  computer  then  asks 

DEVIATIONS? 

If  the  table  of  x,y,  calculated-y , deviation-y  and  RMS 
is  desired,  answer  Y;  otherwise  answer  N 

22.  The  computer  will  ask 

REPLOT? 

If  the  data  is  to  be  modified,  eliminated,  or  extended, 
or  if  another  plot  is  to  be  made  changing  either  x or  y 
limits,  or  the  description,  or  the  fitting  function, 
respond  with  Y;  otherwise  respond  N and  proceed  at  step 
6 

23.  The  computer  will  ask 

NUMBER? 

If  there  are  no  changes  go  to  step  23(c). 

(a)  The  computer  expects  the  index  number  of  the  data 

point  being  changed.  It  will  accept  an  index  number 
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one  greater  than  the  length  of  the  table  (omit  the 
number  corresponding  to  1E35)  with  the  assumption 
that  the  table  is  being  extended.  The  table  may  be 
extended  repeatedly  by  entering  an  index  number  equal 
to  one  more  than  the  previous  maximum  number.  If  the 
number  typed  is  invalid,  step  23  is  repeated.  After 
the  number  is  accepted,  the  computer  will  type 

X,  Y = 

Respond  with  the  value  of  x,  a comma  , the  value  of 
y return.  Step  23  will  be  repeated.  The  data  being 
typed  will  be  logged  according  to  responses  made  in 
steps  7,8 

CAUTION:  DO  NOT  ENTER  1E35  FOR  THE  LAST  NUMBER. 

The  program  knows  how  many  points  there  are  in  the 
table . 

(b)  If  a value  of  x,y  is  to  be  eliminated,  after  its 
index  number  has  been  accepted,  and  X,Y  = has  printed, 
type 

1E35,1E35  return 

(c)  When  all  changes  have  been  made,  answer 

NUMBER? 
with  0 return 

The  procedure  restarts  at  step  14 
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ERRORS 

1.  Alphabetic  errors:  If  an  error  is  discovered  before  re- 
turn is  pressed,  and  the  information  was  alphabetic,  as 
in  steps  7,8,9,16,17,18,19,21,22  type  the  rubout  key; 
the  computer  will  respond  with  a carriage  return;  retype 
the  entire  response,  beginning  with  the  first  character. 

2.  Numeric  errors:  If  an  error  is  discovered  before  a comma 
or  return  is  pressed,  and  the  information  was  a number, 
as  in  steps  10,12,14  and  23,  type  the  rubout  key;  the 
computer  will  not  respond  (automatic  feature  of  the  27 
bit  floating  point  package) ; then  type  the  entire  number 
over  again,  beginning  with  the  first  character. 

CAUTION:  If  two  numbers  were  being  inputted,  separated  by 
a comma,  the  rubout  will  erase  the  CURRENT  NUMBER  ONLY. 

3.  Incorrect  responses  to  Y or  N;  The  computer  will  specifi- 
cally look  for  the  Y.  Any  character  other  than  Y will  be 
treated  as  N except  a carriage  return,  which  should  be 
used  only  after  a character  has  been  typed. 

RESTART  ADDRESS:  0200 
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Figure  29:  PLOT  Sample  Session 
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Fiyure  30:  PLOT  Sample  Plot 


Y=+  2 169642  X*X  - 357142  X + 13  616071 


SAMPLE  DATA 
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7.1  LSFIT 

7.1.1  Program  Description 

Program  LSFIT  was  written  to  calculate  nth  degree  least 
square  curve  fits.  It  will  accept  up  to  1000  pairs  of  x,y 
points  and  computing  the  fitting  function  of  degrees  one 
through  ten.  The  input  is  on  an  a tape  (TAPEl)  and  the  out- 
put is  a listing  of  calculated  values  for  each  degree  and  a 
plot  either  PEN  or  CRT  for  each  degree  fit. 

Note  the  program  is  easily  modified  to  accept  all  data 
from  cards. 
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7.1.2  Control  Cards 

JOBNM,CM65000 , T100 . ID  NO. 

ATTACH , LGO , LSFITBIN , ID=GLASS . 

ATTACH , TAPE1 , DATA, ID=GLASS . 

REQUEST, TAPE 3 9, *Q. 

DISPOSE, TAPE 3 9, *FL . 

ATTACH , CRT , OFFLINECRT . 

LIBRARY, CRT. 

LDSET , PRESET= ZERO . 

LGO. 


NAME 
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7.1.3  Input  Requirements 

RECORD  1 

cc  1-30  execution  time  format  statement  defining  how  the 
x,y  pairs  are  to  be  read. 

ex.  (F10 . 0 , 10X, F12 . 3)  implies  the  x variable  will 
be  found  in  cc  1-10  and  the  y will  be  found  in  cc 
21-32 

RECORD  2 

cc  1-30  data  description  that  will  be  printed  and  plotted 
RECORD  3 OPTION  VARIABLES 

cc  1-15  minimum  x value  (lower  axis  limit) 
cc  16-30  maximum  x value  (highest  axis  limit) 
cc  31-45  minimum  y value  (lower  axis  limit) 
cc  46-60  maximum  y value  (highest  axis  limit) 
cc  61-63  plot  type  PEN  for  pen  plot 

CRT  for  CRT  plot 

cc  66-70  not  used 

cc  71-75  = 1 if  calibration  data 

= 0 if  otherwise 

RECORD  4-N 

DATA  in  the  format  specified  on  record  1 
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7.1.4  Output  Description 

The  first  page  of  output  (figure  31A)  lists  the  input 
x and  y values  in  columns  one  and  two.  The  remaining  ten 
columns  list  the  calculated  least  square  values  for  each 
x value. 

The  second  page  (figure  31B)  lists  the  coefficients  and 
RMS  for  each  fitting  function.  The  first  number  given  is 
the  constant  term,  followed  by  the  slope  and  then  the  quadratic 
coefficient,  etc. 

The  sheet  (figure  31C)  is  a table  of  50  values  equally 
spaced  over  the  x axis  limits  and  the  corresponding  results 
of  applying  each  degree  of  fit. 


A sample  plot  is  shown  in  figure  32. 
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Appendix  PMS-lD  Data  Acquisition  Tape  Description 


Each  record  on  the  PMS-lD  data  acquisition  tape  contains 
four  seconds  of  data.  Each  second  has  data  for  four  probes. 
These  probes  are  subdivided  into  sixteen  24  bit  words.  The 
24  bit  words  are  composed  of  four  6 bit  characters  (four  bit 
BCD  plus  2 fill  bits) . 

The  record  size  (in  bits)  is  then: 


4 seconds  probes  words  bits  = bits 

record  seconds  probe  ~ word  record 


In  CDC  660  words  this  becomes 


bits  1 word 

X - 

record  60  bits 


102.4  words 


This  input  record  format  is  illustrated  in  Appendix  2. 


The  first  of  the  four  probes  contains  VCO  information. 

This  includes  the  PMS  elapsed  second  counter,  the  real  time 
clock  values,  and  the  voltage  levels  for  each  of  the  VCO 
devices.  The  remaining  three  probes  (SCATTER, CLOUD  and  PRECIP) 
each  have  a hosekeeping  word  as  its  first  element.  The  number 
of  particles  counted  in  size  channels  1-15  of  each  device  are 
words  2-16.  The  structure  of  these  four  probes  is  shown  in 
Appendix  3. 

The  housekeeping  word  on  the  probes  (excepting  the  VCO) 
is  a voltage/temp  indicator  to  monitor  key  points  of  the 
particular  device.  As  a means  of  monitoring  ten  different 
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Appendix  1:  PMS-lD  Data  Acquisition  Tape  Description  (cont'd) 

parameters  with  the  same  word,  the  units  of  the  elapsed 
second  clock  determines  which  housekeeping  parameter  is 
being  recorded.  Appendix  4 lists  the  elapsed  second  digit 
and  the  corresponding  housekeeping  reading  for  each  probe. 


Appendix  2:  C130  PMS-lD  Data  Acquisition  Tape  Format 


Record  length:  1024  characters 

Character  length:  4 bit  BCD  + 2 leading  bits 

Buffer  size:  103  60  bit  words  or 

256  24  bit  words 
Time:  4 seconds/record 

64  24  bit  words/second 
Parity:  odd 

Density:  556  bits  per  inch 


6600 

word#  (-*  60  bits  — ►] 


see  Appendix  3 for  de 
scription  of  64  PMS 
(24  bit)  words 


second  1 

6600  words  1-25.6 


second  2 

6600  words  25.6-51.2 
second  3 

6600  words  51.2-76.8 


second  4 

6600  words  76.8-102.4 


elapsed 


See  Appendix  4 for  expansion  of 


App<  .idix  4:  PMS-1D  Status  Word  Definition 
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Appendix  5:  KNOLL ID  TAPE 2 Format 


RECORD  LENGTH: 
TIME: 

PARITY: 

DENSITY: 


130  words  (60  bit) 

1 rec/interval  ave  data 
1 rec/second  full  data 
odd 

variable 


WORD (S) 

VARIABLE 

DESCRIPTION 

1 

TAIR 

A/C  TIME  (HH : MM: SS ) 

2 

TPMS 

PMS  TIME  (HH : MM : SS ) 

at  interval  start  of 

3 

ESEC 

PMS  COUNTS 

ave  data 

4 

ITYPE 

PARTICLE  TYPE 

5 

INTA 

AVERAGING  INTERVAL 

6 

IFORM 

OUTPUT  OPTION 

7 

INTRP 

INTERPOLATION 

8 

FLIGHT  ID 

FLT  E78-01 

9 

FLIGHT  DATE 

DD  MON  YR  (@12@MAR@78) 

@ = blank 

10-24 

DEN (1-15,1) 

DENSITY  (SCATTER) 

25-39 

DEN  (1-15, 2) 

DENSITY  (CLOUD) 

40-54 

DEN (1-15, 3) 

DENSITY  (PRECIP ) 

55-67 

VCO  (13) 

VCO ' s 

68-71 

LWC  ( 4 ) 

LWC  (Sc, Cl, Pr, Total)  / 

see  Appendix  6 for 

72-75 

Z (4) 

Z (Sc, Cl, Pr, Total)  f 

details 

76-79 

DO  (4) 

MEDIAN  VOLUME  DIAMETER 

(Sc, Cl, Pr, Total) 

80-83 

K (4) 

K = M/ (Z) 1/2  (Sc,Cl , Pr 

, Total) 

84 

FF 

FORM  FACTOR 

85 

SFACT 

STABILITY  FACTOR 

86-130 

CTR  (3,15) 

EQUIVALENT  MELTED  DIAMTERS 
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Appendix  6:  KN0LL1D  TAPE2  VCO  Placement 


WORD  KNPLTlD  AXIS  CODE  NAME 


55 

1 

56 

2 

57 

3 

58 

4 

59 

5 

60 

6 

61 

7 

62 

8 

63 

9 

64 

10 

65 

11 

66 

12 

67 

13 

68 

14 

69 

15 

70 

16 

71 

17 

72 

18 

73 

19 

74 

20 

75 

21 

PRESSURE 

EWER 

HEIGHT 

TEMPERATURE 

DEWPOINT 

LWC-JW 

IND  AIRSPEED 

A PRESSURE 

UNUSED 

HEADING 

CALCULATED  AIRSPEED 

TRUE  AIRSPEED 

UNUSED 

LWC  SCATTER 

LWC  CLOUD 

LWC  PRECIP 

LWC  TOTAL 

Z SCATTER 

Z CLOUD 

Z PRECIP 

Z TOTAL 
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Appendix  7:  KNOLL ID  TAPE 3 Format 


TAPE  3 


RECORD 

TIME: 

PARITY: 

LENGTH: 

12  word  (60  bit) 

1 rec/interval 

even  (formatted) 

DENSITY 

variable 

FORMAT  (IX, 

A10, IX, 

4 (1PE10 . 3) , IX 

,4 (1PE10.3) ,1X,  4 ( 1PE10 . 3 ) ) 

WORD (S) 

VARIABLE 

FORMAT 

DESCRIPTION 

1 

ITIM 

A10 

CLOCK 

2-5 

LWC  ( 4 ) 

1PE10. 3 

LWC  ( Sc, Cl, Pr, Total) 

6-9 

Z (4) 

1PE10 . 3 

Z (Sc , Cl, Pr, Total ) 

10-13 

DO  (4) 

1PE10. 3 

DO  (Sc, Cl, Pr, Total) 
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Appendix  8:  KNOLL ID  TAPE 4 Format 


TAPE  4 (PUNCH  FILE) 

RECORD  LENGTH:  8 words  (60  bit) 

TIME:  1 rec/interval 

PARITY:  even  (formatted) 

FORMAT  (1X,3A2,1H. ,5X,F7.1,5X,4 (1PE10.3)) 


WORD (S) 

VARIABLE 

COLUMNS 

DESCRIPTION 

1 

IH 

2-3 

Hours 

2 

IM 

4-5 

Minutes 

3 

IS 

6-7 

Seconds 

4 

HTFT 

14-20 

Height (feet) 

5-6 

LWC (2) 

16-35,36-45 

LWC (PR, Total) 

7-8 

Z (2) 

46-55, 56-65 

Z (PR, Total) 
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Appendix  9:  KNOLLlD  TAPE 7 Format 

TAPE  7 (RAPP  TAPE ) 

RECORD  LENGTH: 

TIME: 

PARITY: 

DENSITY: 

WORD  PARAMETER 

1 DATE  (731222.)  = 22  DEC  73 

2 TIME  ( 24157.)  = 02:41:57 

3 averaging  interval  (seconds) 

4 Water  Content  (gm/m3) (Scatter  Probe) 

5 Water  Content  (gm/m3) (Cloud  Probe) 

6 Water  Content  (gm/m3) (Precip  Probe) 

7 Radar  Reflectivity  (mm6/m3 ) (Scatter  Probe) 

8 Radar  Reflectivity  (mm6/m3) (Cloud  Probe) 

9 Radar  Reflectivity  (mm6/m3 ) (Precip  Probe) 


9 words  (60  bit) 

1 rec/averaging  interval 
odd 

variable 
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Appendix  10:  LEAR  210  Word  Processed  Tape 


RECORD  LENGTH: 
TIME: 

PARITY: 

DENSITY: 


210  words  (60  bit) 
1 record/second 
odd 

variable 


WORD 


PARAMETER 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31-75 

76-120 

121-165 

166-210 


FLAG  fixed  to  7777. 

DATE  (731222.)  = 22  DEC  73 
TIME  ( 24157.)  = 02:41:57 
SEC  fixed  to  1. 

TAS  (m/sec) 

PRES (mb) 

ALT  (km) 

TEMP  ( C) 

Water  content  (gm/m3) (Scatter  Probe) 

Water  content  (gm/m3)  (Cloud  Probe) 

Water  content  (gm/m3) (Precip  Probe) 

Radar  reflectivity  (mm6/m3) (Scatter  Probe) 

Radar  reflectivity  (mm6/m3) (Cloud  Probe) 

Radar  reflectivity  (mm6/m3 ) (Precip  Probe) 

Total  counts  (Scatter  Probe) 

Total  counts  (Cloud  Probe 
Total  counts  (Precip  Probe) 

TWCI  water  content 
TWCI  freql 
Dewpoint  (°C) 

No.  Density  (SZD/m  ) (Scatter  Probe) 

No.  Density  (SZD/m  ) (Cloud  Probe 

No.  Density  (SZD/m  ) (Precip  Probe) 

TWCI  f req2 
TWCI  tempi 
TWCI  temp 2 

largest  particle  size 
average  particle  size 
JW-LWC 

particle  type  code 

counts  for  size  channels  1-45 

unormalized  Number  Density  for  size  classes  1-45 
(SZD/m3) 

water  content  for  size  channels  1-45  (gm/m3) 
equivalent  melted  diamter  for  size  channels  1-45  (u) 


Appendix  11:  LEAR  PMS-lD  Data  Acquisition  Tape  Format 


RECORD  LENGTH: 
CHARACTER  LENGTH: 
BUFFER  SIZE: 

TIME: 

PARITY: 

DENSITY: 


8192  CHARACTERS 

4 BITS  BCD 

547  60  BIT  WORDS  OR 

2048  16  BIT  WORDS 

32  SECONDS/RECORD 

64  16  BIT  WORDS/SECOND 

ODD 

800  BPI  (9  TRACK) 


60  bits 


cloud 


01 


6 


8 


L 


9 


S 


P 


z 


z 


T 


ZZZ6 

X3H 


01 


6 


P 


Z 


Z 


I 


awii 

aasdvaa 


<£> 


aasnNn 

aasnNn 

aasnNn 

aasnNn 

1 

SHflOH 

saNoaas 

SaiDNIW 

p IDMX 

e idmj, 

Z I DM1 

T IDMtL 

3D Vi 70 A 
3Ma-Mr 

aoviioA 

dMaa 

aovxaoA 

awax 

30YI.7OA 

D SS3Hd 

30VI.10A 

ssand 

3618 

X3H 

C 

O 

•r-i 

■P 

<0 

o 

o 


O 

s 

Q 


I 

w 


S 


a, 


TJ 

C 

0) 

CU 

Ci 

< 


CO 

ro  p, 

•H 

o 

0) 

M 

a 


(T> 


O 

U 

> 
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Appendix  13:  PMS-2D  Particle  Image  Tape  Description 


Data  collected  by  the  cloud  and  procip  PMS-2D  data  ac- 
quisition devices  are  written  to  tape  by  a PERTEC  recorder. 

The  data  collected  by  either  device  is  double-buffered  to 
minimize  loss  of  data. 

There  are  two  kinds  of  records  written  to  this  tape, 
slow  and  fast. 

Slow  data  records  are  86  CDC  6600  words  long  and  contain 
the  real  time  clock  value  and  VCO  readings  for  each  second. 

One  slow  record  is  written  every  ten  seconds  and  thus  contains 
ten  seconds  worth  of  data. 

Fast  data  are  the  records  produced  by  the  PMS-2D  buffers, 
and  are  547  60  bit  words  in  length.  These  records  contain 
the  results  of  1024  scans  of  the  32  diodes.  Each  scan  is  a 
binary  representation  of  each  diodes  state.  If  when  scanned 
a diode  was  on  a zero  is  recorded,  if  off  a one  is  recorded. 

In  order  to  minimize  data  written  all  blank  scans  are  counted 
and  when  the  next  particle  is  detected,  the  number  of  blank 
scans  is  written  onto  the  tape.  This  counter  is  prefaced  by 
a bit  pattern  of  10101010  to  differentiate  it  from  a regular 
diode  scan.  The  first  scan  has  the  elapsed  second  counter, 
overload  status  and  a probe  indicator.  If  one  buffer  is  being 
written  to  tape  and  the  other  buffer  for  this  probe  is  also 
filled  the  probe  is  said  to  be  overloaded. 
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Appendix  14:  C130/LEAR  PMS-2D  Fast-Data  Record  Format 


RECORD  LENGTH:  1025  32  bit  scans 

547  60  bit  words 
SCAN  RATE:  1 every  250  nano- 

seconds maximum 


6600  WORD# 


1 

time  & id  word  * 

| scan  1 

2 

| scan  2 

| scan  3 

3 

1 scan 

4 

| scan  5 

4 

| scan  6 

scan  7 

5 

scan  8 

"Tscan  9 

6 

| scan 

o 

7 

scan  11 

j scan  12 

8 

scan  13 

scan  14 

9 

scan  15 

I scan  16 

10 

scan  17 

Iscan  18 

539 

1 scan 

1009 

Iscan  1010 

540 

1 scan  1011 

1 scan  1012 

541 

scan  1013 

542 

scan  1014 

1 scan 

1015 

543 

scan  1016 

I scan  1017 

544 

scan  1018 

scan  1019 

545 

scan  1020 

Iscan  1021 

546 

547 

jscan  1022 
{scan  1024 

scan  1023 

60  bits 


* 


see  appendix  17  for  time  and  id  word  description 


precip  1 2D  clock 
total  ' percent 


Appendix  16:  C130  PMS-2D  Slow  Data  Record  Word  Allocation 

17  18  19  20  21  22  23  24  25  26  27  28  29  30  31  32  word# 
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Appendix  17:  C130/LEAR  Time  & ID  Word  Description 


TIME  BYTES 


byte  # 

1 
2 

3 

4 

— 8 bits  

note:  every  "seconds"  digit  is  4 bit  BCD 


elapsed  seconds  x 1000 

elapsed  seconds  x 100 

elapsed  seconds  x 10 

elapsed  seconds  x 1 

elapsed  seconds  x 0.1 

elapsed  seconds  x 0.01 

elapsed  seconds  x 0.001 

Id  4 Id  3 Id  2 Id  1 

IDENTIFICATION  BITS 


FAST  DATA 

SLOW 

C130E 

C130A 

DATA 

Id 

4 

data 

rH 

source 

O 

O 

in 

Id 

3 

data 

o 

<U 

C 

source 

X 

o 

Id 

2 

overload 

w 

a 

o 

rH 

rH 

u 

Id 

1 

overload 

<D 

(/) 

Id  4 = I# record  from  cloud  size  probe 
Id  3 = If  record  from  precip  size  probe 
Id  2 = 1,  cloud  size  probe  overload 
Id  1 = If  precip  size  probe  overload 

note:  Id  4 and  Id  3 may  not  be  one  in  the  same 
record 


265 


Appendix  18:  PMS-2D  Particle  Tape  Format 


RECORD  FORMAT:  SCOPE-NOS/BF 

RECORD  LENGTH:  VARIABLE  MINIMUM  12  WORDS,  MAXIMUM  512 
RECORD  TYPE:  BINARY 


WORD  DESCRIPTION 

1 
2 

3 

4 

5 

6 

7 

8 
9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 


512  VOLUME  IN  CUBIC  DIODES 


Each  record  consists  of  the  12  word  identification  block  and 
then  up  to  50  particle  descriptions,  each  10  words  long.  One 
fast  input  record  to  TWODEE  can  generate  one  or  more  of  these 
records.  The  length  of  each  output  record  is  determined  by 
the  number  of  particles  in  each  fast  input  record. 


TIME  OF  RECORD  IN  FORM  HH:MM:SS.F  (DISPLAY  CODE) 
SLOW  RECORD  AND  FAST  RECORD  NUMBER  IN  FORM 
WORD  2 = SLOW* 10 000  + FAST  (BINARY) 

SAMPLE  TIME  IN  MS  (BINARY) 

CLOCK  SAMPLING  RATE-PERCENTAGE  (BINARY) 

PROBE  (BINARY) 

OVERLOAD  INDICATOR  (BINARY) 

KISTLER  AND  BACKUP  PRESSURE  IN  FORM 
WORD  7 = KIST*10000  + PRESS  (BINARY-VCO  COUNTS) 
PRESSURE  GRADIENT  VCO  COUNTS (BINARY-VCO  COUNTS) 
TEMPERATURE  VCO  COUNTS  ( " ) 

DEWPOINT  VCO  COUNTS  ( " ) 

TRUE  AIRSPEED  ( " ) 

JW-LWC  VCO  COUNTS  ( " ) 

AREA  IN  SQUARE  DIODES -NEGATIVE  IF  EDGE  REJECTION 

(BINARY) 

PERIMETER  IN  DIODES  ( " ) 

HORIZONTAL  FERET  PROJECTION  IN  DIODES  (BINARY) 


VERTICAL  FERET  PROJECTION  IN  DIODES  ( " ) 
HORIZONTAL  PROJECTIONS  IN  DIODES  ( " ) 
VERTICAL  PROJECTIONS  IN  DIODES  ( " ) 
UNUSED 

MAXIMUM  LENGTH  IN  DIODES  ( " ) 
THETA  ANGLE  OF  PARTICLE  ORIENTATION  ( " ) 
VOLUME  IN  CUBIC  DIODES  ( " ) 
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Appendix  19:  SPANDAR  Input  Tape  Format 


TAPE1 

PARITY: 

DENSITY: 

RECORD: 

TIME: 


(SPANDAR) 

EVEN 

800  BPI  (7  TRACK) 
VARIABLE  LENGTH 
1 SECOND  PER  DATA  RECORD 


The  SPANDAR  input  tape  consists  of  multiple  sets  of: 
1 header  record  followed  by  many  data  records 


HEADER  RECORD  (28  bytes) 


BYTES  TYPE 


DESCRIPTION 


6 

2 

2 

2 

2 

2 

2 

2 

4 

4 


8 bit  EBCDIC 
fixed  binary 

It 


floating  point 

II 


A/C  type  (6  char) 

Run  ID  Number  (pass) 

Day  of  Scan  (DOY) 

Year  of  Scan  (2  digits) 
Scan  Number 
Time  of  Scan  (hours) 
Time  of  Scan  (minutes 
Number  of  data  records 
Time  of  Scan  (seconds) 
Nominal  Altitude  (feet) 


DATA  RECORDS  one  per  second  (40  bytes) 


2 fixed  binary 

2 

4 floating  binary 

4 
4 
4 
4 
4 
4 
4 
4 


Greenwich  Time  (hours) 
Greenwich  Time  (minutes) 
Greenwich  Time  (seconds) 
Average  Reflectivity  (dbz) 
Elevation  (degrees) 

Azimuth  (degrees) 

Radar  Slant  Range  (nm) 
Ground  Range  (nm) 

Ground  Range  (nm) 

Altitude  (feet) 

Altitude  (km) 
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Appendix  20:  SPANDAR  TAPE 3 Format 


TAPE  3 

TYPE:  SCOPE-NOS/BE  STANDARD 

RECORD:  VARIABLE  LENGTH 

TIME:  1 SECOND  PER  DATA  RECORD 

HEADER  RECORD  LENGTH  2 WORDS 

WORD  1 ZERO  OR  DAY  OF  YEAR 

WORD  2 FOR  SPANDAR  TAPE 

DATA  RECORD  LENGTH  7 WORDS  (DISPLAY  CODE) 

SPANDAR  output  tape  consists  of  multiple  sets  of  header 
record  followed  by  a variable  number  of  data  records 


FORMAT  DESCRIPTION 


2X 
13 
13 
F7.3 
F9 .0 
F8 . 3 
F7 . 3 
F9.0 
F7.1 
F6 . 0 


blank 

GMT  (hours) 

GMT  (minutes) 

GMT  (seconds) 

radar  slant  range  (meters) 
azimuth  (degrees) 
elevation  (dgrees) 
altitude  (meters) 
reflectivity  (dbz) 
counts  (forced  to  zero) 
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Appendix  21:  SPANDAR  TAPE4  (Punch  File)  Format 

The  SPANDAR  output  deck  consists  of  multiple  sets  of: 
1 header  card  followed  by  many  data  cards. 

HEADER  CARD 

DESCRIPTION 

1 

Date 

Aircraft 
Pass  No. 

Altitude  (feet) 

Scan  No. 

No.  of  Data  cards  that  follow 
Start  Time  (hhmmss.f) 

Stop  Time  (hhmmss.f) 

DATA  CARDS 

DESCRIPTION 

Blank 

Time  (hhmmss.f) 

Blank 
Z (DBZ) 

Elevation  (deg) 

Azimuth  (deg) 

Radar  Slant  Range  (NM) 

Ground  Range  (NM) 

Ground  Range  (NM) 

Altitude  (feet) 

Altitude  (km) 


COLUMNS 

1 

2-9 

10-14 

15-22 

23-30 

31-38 

39-46 

47-54 

55-62 

63-70 

71-78 


COLUMNS 

1 

4-13 

16-21 

24-26 

29-36 

39-41 

44-46 

49-56 

59-66 
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Appendix  22:  MORT  Input  Tape  Format 


HEADER  RECORD 


5A6 

Title 

E15 . 7 

RCS  o-count  equivalent  (dBsm) 

E 1 5 . 7 

RCS  count  increment  equivalent 

(dBsm) 

A6 

Radar  sensor  used 

RECORD 

213, F7. 3 

Time  (Hours,  minutes,  seconds) 

F9.0 

Range  (meters) 

F8 . 3 

Azimuth  (degrees) 

F7. 3 

Elevation  (degrees) 

F9 . 0 

Altitude  (degrees) 

F7 . 1 

Z (dB-mm6 /m3 ) 

F6. 1 

RCS  counts 

F8.3 

M (gm/m3 ) 

F6 . 1 

Polarization  Ratio  (PP  to  OP  in 

dB) 

F6 . 1 

Average  S/N  (dB) 

F6.1 

Average  cross  section  (dBsm)  = 

. 55118 

15 

Total  number  of  data  points  available 
averaging  interval 

15 

Number  of  points  accepted  for  use 

counts 

in 


170 
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Appendix  23 

VARIABLE 

HR :MN : SC . F 

Z (DBZ) 

EL  (DEG) 

ZA  (DEG) 

RSLRA (NAM) 
GRRA(NM) 

GRRA (KM) 

ALT (FT) 

ALT (KM) 


: RAPP  Radar  Tape  Format 

DESCRIPTION 

HOURS,  MINUTES  AND  DECIMAL  SECONDS 

RADAR  READING 

ELEVATION  OF  RADAR  DEGREES 

AZIMUTH  OF  RADAR  DEGREES 

SLANT  RANGE  NAUTICAL  MILES 

GROUND  RANGE  NAUTICAL  MILES 

GROUND  RANGE  KILOMETERS 

ALTITUDE  IN  FEET 

ALTITUDE  IN  KILOMETERS 


VCO  CALIBRATION 

CHANNEL  SYMBOL  DESCRIPTION  INTERCEPT  SLOPE  3RD  TERM  UNITS 
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Appendix  24:  C130E  VCO's  (PMS-lD  system) 


vco 

CHANNEL  SYMBOL  DESCRIPTION  CALIBRATION  UNITS 
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